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FOREWORD 


This  report  (Litton  Publication  No.  7244)  was  prepared  by  the  Space 
Sciences  Laboratories  of  Litton  Systems,  Inc.  ,  at  336  N,  Foothill  Road, 
Beverly  Hills,  California,  under  Contract  No.  AF  33(6 15) -3744,  cover¬ 
ing  the  period  from  1  July  1966  to  31  July  1967. 

This  document  is  a  revision  of  prior  Litton  publication  No.  6730, 
and  was  monitored  by  the  Air  Force  Aero  Propulsion  Laboratory  with 
Mr.  Richard  B.  Rivir  as  the  USAF  Project  Engineer. 

Publication  of  this  report  does  not  constitute  Air  Force  approval  of 
the  report's  findings  or  conclusions.  It  is  published  only  for  the  exchange 
and  stimulation  of  ideas. 


Paul  E.  Peko,  Major,  USAF 
Chief,  Propulsion  and  Power  Branch 
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ABSTRACT 


The  work  described  in  this  report  comprises: 

a.  Theoretical  studies  of  the  performance  of  an  electromagnetic 
induction  plasma  accelerator  for  radiation-cooled  applications, 
and  for  pulsed  operation. 

b.  The  improvement  of  pre-existing  instrumentation  for  the  measure¬ 
ment  of  developed  thrust  and  input  electric  power. 

c.  The  development  of  suitable  materials  and  techniques  for  the 
construction  of  thrustors. 

d.  The  experimental  testing  of  thrustors  for  the  purpose  of  measuring 
performance  and  life  times.  ' 
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INTRODUCTION 


The  work  described  in  this  report  i6  a  continuation  of  work  done 
under  Contract  AF  33(657)8697  and  summarized  in  Report  AFAPL-TR- 
65-491,  which  will  be  referred  to  hereinafter  as  the  Report,  since  many 
occasions  will  be  found  to  use  information  therein  presented. 

The  work  here  described  can  be  divided  in  four  major  areas: 

a.  Theoretical  studies  of  the  performance  of  the  thrustor  for  the 
case  in  which  ferrites  cannot  be  used  because  the  operating 
temperature  exceeds  the  curie  point  of  the  materials.  This 
condition  applies  when  the  thrustor  is  operated  in  the  radiation- 
cooled  mode.  Another  theoretical  aspect  investigated  was  pulsed 
operation,  in  which  high  instantaneous  power,  compatible  with 
high  efficiency,  is  used  for  short  periods  of  time  so  as  to  generate 
low  average  thrust  at  low  average  power  level. 

b.  Improvement  of  the  experimental  facility,  particularly  in  respect 
to  the  maintaining  of  better  vacuum  during  long  thrustor  firing, 
and  the  improvement  of  the  precision  of  measurement  for  thrust 
and  electrical  input  power. 

c.  Development  of  materials  and  techniques  for  the  fabrication  of 
thrustors  capable  of  extended  firing  times  in  the  radiation-cooled 
mode  as  well  as  in  the  water-cooled  mode. 

d.  Experimental  measurement  of  performance  data,  conditioned 
upon  the  availability  of  thrustors  capable  of  being  operated  for 
sufficiently  long  times  at  sufficiently  high  *:uwer  levels,  and 
therefore  upon  satisfactory  progress  under  point  c.  above. 

Work  done  in  each  of  the  above  areas  is  described  hi  detail  in  the 
remainder  of  this  report. 


1  Fonda- Bonardi,  G.  "Development  of  a  plasma  accelerator  for  thrust,  " 
Final  report  on  Contract  AF  33(657)8697,  Litton  Publication  No.  42  31, 
April  1 967 
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THEORETICAL  CONSIDERATIONS 


2.  1  GENERAL 


Three  topics  have  been  considered  during  the  present  research 
period.  The  first  two  concern  the  present  experimental  device  and  the 
third  topic  considers  a  possible  future  ihrustor  mode  of  operation.  In 
the  orde  <  scussed  they  are: 


2.  2  Radiation-cooled  thrustor  design  and  performance 

2.  3  Thrustor  models 

2.  4  Pulsed-mode  thrustor  operation 

One  of  the  main  objectives  of  this  research  period  was  the  develop¬ 
ment  of  a  radiation-cooled  thrustor  in  contrast  to  the  water-cooled  thrus¬ 
tor  of  Contract  AF  33(657)-8697.  Experiments  showed  (see  Section  5.2) 
that  at  the  design  power  levels,  a  radiation  cooled  coil  would  run  at 
temperatures  exceeding  the  Curie  point  of  magnetic  materials  and  far 
beyond  that  of  ferrites.  A  radiation -cooled  coil  was  designed  by  methods 
shown  in  the  Report  sec  2  and  details  are  given  below.  The  coil  has 
closely  the  same  dimensional  and  magnetic  characteristics  of  the  pre¬ 
vious  water-cooled  coils  at  the  expense  of  higher  I^R  losses.  The  pre¬ 
sent  performance  analysis  includes  these  losses.  Since  the  original 
work  was  done  at  240  Khz  and  all  the  instrumentation  was  developed  for 
this  frequency,  a  performance  vs.  frequency  analysis  was  made  to  see 
whether  performance  peaks  were  close  enough  to  240  Khz  or  whether  a 
large  shift  might  be  needed. 


The  losses  associated  with  fractional  n 


analyzed  for  a  thrustor  model  in  which  resonant  charge-exchange  effects 
are  not  dominant,  and  an  estimate  was  made  of  the  charge-exchange 
mean  free  path  in  the  interaction  region  of  the  thrustor. 
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The  high  efficiencies  (indicated  in  the  Report)  of  the  water-cooled 
thrustor  were  obtained  at  high  power  levels,  which  may  exceed  the  average 
power  available  for  certain  missions  for  which  the  thrustor  may  be  other¬ 
wise  attractive.  However,  if  the  thrustor  could  be  run  in  a  pulsed  mode 
at  low  duty  cycle,  high  peak  power  could  be  applied  at  low  average  power, 
therefore  taking  advantage  of  the  high  power-high  efficiency  thrustor 
characteristics.  At  the  time  of  this  writing,  high  efficiency  data  points 
were  not  available  for  the  current  coil  design  so  the  original  data  points 
of  the  Report  were  used. 

2„  2  RADIATION- COO  LED  THRUSTOR  DESIGN  AND  PERFORMANCE 
CALCULATIONS 

Figures  2-1,  2-2,  and  2-3  show  respectively  the  magnetic  flux  pattern, 
the  inductive  potential  pattern,  and  the  axial  dependence  of  the  inductive 
potential  for  the  radiation-cooled  thrustor.  These  are  to  be  compared 
to  F;<mres  2-23,  2-24  and  2-13  of  the  Report  which  are  the  corresponding 
plots  for  the  water-cooled  thrustor.  It  should  be  noted  that  the  Report 
calculation  were  done  for  coil  currents  in  ampereturns  and  the  present 
calculation  expresses  the  current  in  amperes.  The  trajectory  equation  for 
the  plasma  (Eqn.  2.  11  of  the  Report)  is: 


,2  7  t  2  -34.  I2z  .  2 

d  z  =  4.355  X  10  I  e  sin  <ut 


where  lg  is  the  coil  current  in  peak  amperes.  The  digital  computer 
available  in  the  laboratory  was  programmed  to  solve  this  differential 
equation  for  a  series  of  value  of  the  operating  frequency  f  =  u/2jt.  The 
first  plot  of  this  series  is  shown  in  Figure  2-4,  which  presents  the  ter¬ 
minal  exhaust  velocity  of  the  thrustor  as  a  funct;on  of  frequency  and 
coil  current.  This  plot  is  the  limiting  case  of  zero  mass  flow  (therefore 
zero  thrust),  i.  e.  there  are  no  de-coupling  effects  due  to  plasma  currents 
and  in  a  way  this  represents  an  ideal  thrustor.  When  plasma  currents 
de-couple  the  coil  from  the  plasma  more  coil  current  is  needed  to  pro¬ 
duce  the  same  output  velocity,  Section  2.  3.  2  of  the  Report  is  an.  analysis 
of  the  loading  effect.  Equation  2.  37  of  the  Report  allows  one  to  calculate 
the  currents  needed  when  the  plasma  loads  the  thrustor  and  Figures  2-5 
through  2-15  show  the  results  of  the  calculation  applied  to  the  radiation- 
cooled  coil.  This  series  show  lines  of  constant  coil  current  as  a  function 
of  frequency  and  exhaust  velocity  for  a  series  of  values  of  thrust  from 
0.  05  to  5.  0  newton. 
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Figure  2-2.  Equipotentials  of  Radiation  Cooled  Thrustor 
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Figure  2-9.  Current  vs  Terminal  Velocity  and  Frequency  {T  =  0.  5  newton) 
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Figure  2-10.  Current  vs  Terminal  Velocity  and  Frequency  {T  =  0.7  newton) 
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Figure  2-11.  Current  vs  Terminal  Velocity  and  Frequency  (T  =  1 . 0  newton) 
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The  power  dissipated  in  the  coil  by  the  I  R  losses  can  be  computed 
if  the  value  of  the  coil  resistance  is  known  as  a  function  of  frequency. 

This  was  comouted  by  the  standard  skin- resistance  formula  in  which  the 
resistivity  of  the  copper  was  taken  to  be  2.  5  times  that  at  room  tempera¬ 
ture  in  accordance  with  the  experimental  results  of  section  5.2.  The 
total  resistance  of  a  coil  having  the  dimensions  shown  in  the  analog  plots 
of  Figures  2-1  and  2-2,  built  of  copper  tube  1/8"  in  diameter  and  operated 
at  550°C  is  given  by 


R  =  0.  229 sfT 

where  f  is  the  frequency  of  operation  in  Megahertz.  The  power  dissipated 
is  therefore 


W  =  -5—1  2  0.  11457TI  2 

2  c  v  o 

where  IQ  is  the  peak  coil  current  in  amperes.'  Using  this  and  the  per¬ 
formance  curves  of  Figure  2-4  through  2-15  a  new  set  of  curves  was 
derived,  showing  power  loss  in  the  coil  as  a  function  of  frequency,  ex¬ 
haust  velocity  and  thrust.  These  are  presented  in  Figure  2-16  through 
2-27.  The  series  of  figures  from  Figure  2-28  through  2-35  is  the  same 
information  in  which  lines  of  constant  power  are  plotted  in  a  thrust- 
velocity  plane  for  a  series  of  frequencies  ranging  from  30  Kilohertz  to 
1  megahertz.  These  graphs  illustrate  the  predicted  performance  of  the 
thrustor  at  constant  frequency.  From  these  it  is  possible  to  compare 
the  dissipated  power,  Wr,  with  the  beam  power  =  Tv/2  where  T  is 
the  thrust  and  v  the  exhaust  velocity.  We  define  an  ohmic  efficiency  r) 
by 


11  r  =  W  +  W,  ’ 
r  k 

The  series  of  figures  from  Figure  2-36  through  2-43  show  lines  of  con¬ 
stant  »ir  in  the  thrust-velocity  plane  for  a  series  of  frequencies.  There 
are  several  points  to  be  noted  concerning  this  last  set  of  curves. 

a.  The  efficiency  q  is  calculated  on  the  basis  of  a  coil  assumed  to 
be  operating  at  a  constant  temperature  of  550°C.  Because  of  the 
dependence  of  resisitivity  on  temperature  the  coil  will  be  cooler 
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Figure  2-16.  Ohmic  Power  Loss  in  Coil  vs  Terminal  Velocity  and  Frequency 

(Zero  Loading:  T  =  0) 
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Figure  2-17.  Ohmic  Power  Loss  in  Coil  vs  Terminal  Velocity  and  Frequency 

(T  =  0.  05  newton) 
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Figure  2-21.  Ohmic  Power  Loss  in  Coil  vs  Terminal  Velocity  and  Frequency 

(T  =  0.  5  newioti) 
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f  igure  ?-22.  Ohmic  Power  Loss  in  Coil  vs  Terminal  Velocity  and  Frequency 

(T  =  0.  7  newton) 
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Figure  2-2  j.  Ohmic  Puwer  Lugs  in  Cuil  vs  Teiiiiiuui  Velocity  and  r  requency 

(T  =  1.0  newton) 
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Figure  2-24.  Ohmic  Power  Loss  in  Coil  vs  Terminal  Velocity  and  Frequency 

(*T  a  1.5  newton) 
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Figure  2-26.  Ohmic  Power  Loss  in  Coil  vs  Terminal  Velocity  and  Frequency 

(T  =  3.  0  newton) 
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Figure  2-29.  Ohmic  Power  Loss  in  Coil  vs  Thrust  and  Terminal  Velocity  at 

Constant  Frequency  (£  =  50  Kilohertz! 
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Figure  2-31.  Ohmic  Powe 
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Figure  2-33.  Ohmic  Power  Less  in  Coil  vs  Thrust  and  Terminal  Velocity  at 

Constant  Frequency  {f  =  500  Kilohertz) 
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Figure  £-54.  Ohmic  Fower  juoss  in  Coil  vs  ihrusi  and  lerminai  Velocity  at 

Constant  Frequency  (i  -  700  Kilohertz) 
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Figure  2-37.  Ohmic  Efficiency  of  Coil  vs  Thrust  and  Terminal  Velocity  at 

Constant  Frequency  {i  =  50  Kilohertz) 
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Figure  2-38.  Ohmic  Efficiency  of  Coil  vs  Thrust  and  Terminal  Velocity  at 

Constant  Frequency  (f  =  70  Kilohertz) 
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Figure  2-41.  Ohmic  Efficiency  of  Coil  vs  Thrust  and  Terminal  Velocity  at 

Constant  Frequency  (f  =  500  Kilohertz) 
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and  therefore  more  efficient  at  lower  dissipation  levels,  say  < 

1000  watts.  For  power  £1000  watts  the  temperature  will  be 
higher,  and  the  losses  relatively  more  severe  which  would  pro¬ 
bably  result  in  thermal  coil  damage.  For  this  reason  the  curves 
are  shown  dotted  for  £  1000  watt. 

b.  The  only  losses  considered  in  the  definition  of  hr  are  the  I^R 
losses  in  the  coil;  *lr  does  not  include  ionization  losses  or  the 
mass  utilization  efficiency.  Therefore  it  is  not  a  suitable  over¬ 
all  efficiency,  but  it  was  computed  for  the  radiation-cooled  coil 
since  the  I^R  losses  are  relatively  more  important  than  in  the 
water-cooled  coil  and  these  losses  should  be  minimized  regard¬ 
less  of  the  presence  of  other  loss  mechanizms. 

c.  From  the  performance  curves  (Figure  2-5  through  2-15)  it  can  be 
seen  that  performance  generally  peaks  between  100  Khz  and  1  Mhz. 
Since  the  instrumentation  developed  under  the  previous  contract 
was  designed  for  240  Khz  it  was  decided  that  the  frequency  should 
remain  the  same.  It  is  not  a  trivial  job  to  change  instrumentation 
and  re-calibrate  for  a  different  frequency;  this  would  only  be 
done  if  results  significantly  better  would  be  expected. 

Figure  2-44  is  a  performance  curve  at  240  Khz  in  the  T -v  plane 
where  the  lines  of  constant  coil  current  are  g;ven  in  rms  amperes. 

2.3  THRUSTOR  MODELS 

The  performance  computed  in  Section  2.2  is  based  on  the  assumption 
that  the  thrust  is  generated  by  the  acceleration  of  ionized  mass  under  the 
action  of  an  electromagnetic  field.  The  current  needed  to  generate  the 
field  is  known,  and  the  I^R  losses  can  be  computed  for  a  given  coil  geo¬ 
metry.  To  these  losses  one  must  add  the  amount  of  power  needed  to 
ionize  the  gas.  The  ionization  energy  of  mercury  is  10.4  electron 
volts,  but  in  general  more  electron  volts  arr  required  to  generate  an 
ion  pair;  typically  several  hundred  electron  volts.  An  efficient  ioniza¬ 
tion  mechanism  must  ionize  a  large  fraction  of  the  number  of  atoms 
present,  and  spend  lew  electron-volts  per  ion  to  do  it.  Une  cannot  pre¬ 
dict  accurately  in  advance  what  fraction  will  be  ionized  in  the  thrustor, 
nor  how  much  power  will  be  used  in  doing  it;  it  is  possible  however  to 
arrive  at  some  estimates  of  the  effectB  on  the  predicted  efficiency  of 
either  loss  mechanism. 


T  (newton) 


50 


. 


OT  30-4* 


If  the  effluent  mass  is  only  partially  ionized,  and  the  ionised  fraction 
is  accelerated  to  the  predicted  velocity,  the  neutral  fraction  may  be  par¬ 
tially  dragged  along  by  collisions  and  charge-exchange  processes.  The 
effect  will  be  strong  if  the  mean  free  path  is  short  compared  with  the 
dimensions  of  the  accelerating  region,  weak  otherwise.  In  the  limit, 
when  the  mean  free  path  is  very  large,  the  effect  will  be  negligible  and 
the  entire  thrust  will  be  due  to  the  ionized  fraction  accelerated  to  the 
full  exhaust  velocity,  while  the  neutral  fraction  will  contribute  nothing. 

In  this  limiting  case,  since  the  effective  velocity  is  computed  dividing 
the  measured  thru6t  by  the  total  mass  flow,  the  observed  effective  vel¬ 
ocity  v  will  be  decreased  by  the  ratio  of  ionized  mass  to  the  total  mass, 
i.  e. ,  by  the  mass  utilization  qm,  and  similarly  for  the  power  efficiency 

v  =  T/m  =  t!mv^ 

l 

i 

1  =  T2/2m  W  =  t]  n 

e  m  * 

n  -  m./m  =  m./m 

m  *  i 

where  the  asterisk  refers  to  the  ionized  portion  only.  This  model  can  be 
called  the  zero-coupling  model,  and  can  be  expected  to  be  valid  at  a  very 
low  value  of  mass  flow,  when  the  density  in  the  interaction  region  is  low. 
It  is  therefore  important  to  compute  the  particle  density  in  the  interaction 
region  as  a  function  of  mass  flow.  Figure  2-45  shows  the  computed  den¬ 
sity  (particles  per  cubic  meter)  as  a  function  of  axial  distance,  z  (cm) 
from  the  nozzle,  under  the  following  assumptions: 

a.  A  line  source  at  the  location  of  the  nozzle 

b.  Exit  velocity  equal  to  sonic  (thermal)  velocity  at  the  expected 
temperature  of  the  nozzle  (260  m/s) 

c.  Cosine  distribution  with  angle  from  the  axis 

Superimposed  on  the  graph  is  a  duplication  of  Figure  2-3,  redrawn  to  the 
new  axial  scale  z,  showing  the  inductive  potential  and  indicating  the  ex¬ 
tent  of  the  region  over  which  the  acceleration  forces  arc  effective.  Also 
shown  are  (dashed)  lines  relating  particle  density  to  charge -exchange 
mean  free  path  for  three  values  of  ion  velocity2;  mean  free  path  lengths 
are  read  in  centimeters  on  the  abscissa.  It  can  be  seen  that  the  mean 

2Hasted,  J.  B.  "Physics  of  Atomic  Collisions",  Washington  1964,  p.  416 
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Figure  2-45.  Particle  Density  in  Interaction  Region 


r  =  z  (cm) 
\  (cm) 


free  path  i?  of  the  order  of  the  dimensions  of  the  interaction  region  for 
all  values  of  mass  flow  of  interest,  so  that  the  zero-coupling  model  may 
provide  a  reasonably  good  approximation  for  the  operation  of  the  thrustor. 
To  be  sure,  collisions  and  charge  exchange  processes  will  occur  but  will 
not  be  dominant;  and  comparison  with  the  experimental  results  will  show 
their  relative  importance. 

A  consequence  of  this  situation  is  chat  the  operating  efficiency  is 
directly  proportional  to  the  fraction  of  ions  produced,  and  any  mechanism 
that  can  increase  it  will  correspondingly  improve  the  efficiency.  There 
are  three  main  mechanisms  for  ion  production; 

a.  The  induction  field,  where  the  loop  voltage  V  is  equal  to  the  flix 
rate  of  change  d$/dt  and  the  electric  field  strength  Eg  is  equal 
'■'•o  V/2ffR.  From  Figure  2-1  we  have,  ai  240  Khz  and  1  cm  from 
the  plane  of  the  coil, 

Eg  =  8.8  I 

where  Eg  ie  in  volt/meter  and  1^  is  ampere. 

b.  The  field  due  to  the  end-to-end  voltage  of  the  coil.  This  field 
may  be  locally  stronger  than  Eg  and  may  provide  some  of  the 
ionization,  but  it  also  provides  a  mechanism  for  ion  removal 
since  it  can  drive  electrons  into  the  walls  of  the  coil.  Since 
these  electrons  then  tend  to  neutralize  the  field,  and  since  the 
maximum  path  length  available  to  this  field  is  short  for  the 
prevailing  density  (the  conditions  are  far  up  the  low-pressure 
branch  of  Paschen's  curve),  the  contribution  of  this  field  is 
probably  minor. 

c.  The  pre-ionization  of  the  effluent  gas  from  the  nozzle.  This 
derives  from  back- streaming  electrons  and  is  always  present 
during  normal  operation  of  the  thrustor.  The  gas  coming  out 
of  the  nozzle  is  strongly  luminous  and  appears  ionized  on  the 
spectrograph.  The  contribution  of  this  source  amounts  probably 
to  a  few  percent. 

The  only  mechanism  capable  of  producing  massive  ionization  ie  the 
first  one.  Sudden  and  large  improvement  of  performance,  observed  in 
other  electrodeless  induction  devices,  were  traced  to  the  dominance  of 
this  mechanism  when  the  Eg  field  becomes  sufficiently  strong  and  causes 
an  avalanche  increase  of  ionization  around  the  loop  of  plasma.  It  is 
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therefore  important  to  estimate  the  range  of  parameters  over  which 
avalanche  ionization  can  take  place. 

Under  steady  state  operation  the  rate  of  production  of  ions  mast  be 
equal  to  the  rate  of  removal,  and  both  must  be  uniform  around  the  plasma 
ring;  the  electron  density  (equal  to  the  ion  density)  must  be  sufficient  to 
insure  that  more  electrons  are  made  und-  r  the  influence  of  the  Enfield. 
The  number  of  electrons  that  one  electron  makes  per  unit  length  is  pro¬ 
portional  to  the  number  of  neutral  atoms  available  (and  therefore  to  the 
density  and,  in  the  case  of  the  thrustor,  to  the  neutral  mass  flow).  In 
the  usual  case  this  is  written  as 


-3 -  =  an  (n  -  n  ) 

dx  e  e 

where  ais  known  as  Townsend's  first  coefficient  and  n  is  the  initial  atom 
density.  In  the  case  of  the  thrustor  the  rate  of  electron  production  is 
equal  to  the  rate  of  ion  ejection  flow  per  unit  length.  By  indicating  time 
rates  with  the  usual  dot  we  have 


— -  =  ah.  (n  -  n  ). 

dx  i  » 

Now  the  rote  of  ion  ejection  per  unit  length  of  plasma  ring,  being  uniform, 
is  equal  to  the  total  rate  of  ion  ejection  di'  led  by  the  total  length  0/  the 
ring;  multiplying  by  the  mass  of  the  ions  have 

rh. 

w  ‘  ori,i (n  -  "i1' 

to  this  we  may  add  a  sma],l  term  b  which  represents  the  contribution  of 
other  sources.  We  have  already  noted  that  the  density  n  is  proportional 
to  th  mass  flow,  as  illustrated  in  Figure  2-45,  say 

n  =  c  m,  n,  =  c  tv 
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we  have  then 


2ffR 


c  a  m  ( m  -  m . )  +  b 
1  \ 


’m  .  ,,  _  .  ,  b 

~r=rr  =  c  an  m  (1  -  r\  )  +  -r- 
2JTK  m  m  m 


2 

n  + 

m 


1 


2jtR  cam 


-  1 


m  .  2 

cam 


■■=  0 


If  no  other  sources  of  electrons  are  present,  so  that  b  =  0,  t|  will  be 
larger  than  zero  u 

2ffRcam>l , 

otherwise  the  plasm*  will  not  sustain  ionization  under  the  action  of  the 
Eg  field.  If  b  is  not  zero  there  will  be  some  ionization  for  lower  values 
of  the  E$  field,  but  the  transition  to  high  density  will  occur  when  the  Eg 
avalanche  takes  over.  Since  a  is  a  function  of  E  and  n  (or  c  m),  taking 
the  published  values3of  a  it  is  possible  to  compute  the  boundary  for 

2*rR  cam  =  I 

as  a  function  ->f  1^  and  m.  This  is  shown  in  Figure  2-46,  The  curve  is 
probably  mor  nearly  correct  for  m<0.  1  grams/sec,  where  the  mean 
free  path  approaches  the  dimensions  of  the  interaction  region  and  the 
zero-couplinc.  model  is  more  nearly  verified.  The  location  of  the  ava¬ 
lanche  thresl. old  becomes  more  uncertain  for  larger  mass  flows,  but  the 
general  shape  of  the  curve  is  probably  correct.  For  comparison,  the 
dashed  line  on  the  same  graph  shows  the  location  of  the  threshold  for  the 
thrustor  with  ferrites  used  in  Contract  AF  33(657)-8697,  subject  to  the 
same  type  of  uncertainty.  The  value  of  b  cannot  be  predicted  a  priori 
except  to  the  extent  that  it  is  likely  to  be  small,  so  as  to  yield  values  of 
^m  between  G.  05  and  0,  10  in  the  region  below  the  avalanche  boundary  to 
conform  with  the  experimental  evidence. 


3  Brown,  S,  C. ,  "Basic  data  of  Plasma  Physics"  Wiley,  New  York  1959, 
p.  129 
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igure  2-46,  Probable  Boundary  of  High  Performance  Area 


Figure  2-47  is  a  plot  of  thrust/current  vs  m  derived  from  the  theo¬ 
retical  curves  for  the  thrustor  at  240  khz.  The  curve  labeled  q  =  1  is 
the  curve  directly  derived  from  the  theory.  The  lines  fanning  out  from 
the  band  are  labelled  by  t.  e  thrust;  these  variations  are  due  to  the  loading 
factor,  j.  e.  the  henz's  law  effect  of  the  plasma  currents  on  the  field  pro¬ 
ducing  them.  The  other  curves  {abbreviated  to  the  range  of  the  present 
data)  are  derived  by  assuming  that  the  same  thrust  is  generated  by  the 
same  number  of  particles  accelerated  for  different  values  of  added  neu¬ 
tral  mass  flow,  so  they  are  labelled  by  mass  utilization  efficiencies.  The 
graph  is  reproduced  in  Section  6.  2  with  the  experimental  points  added, 
where  the  r;m  are  interpreted  as  degrees  of  ionization  and  the  zcrc-coupling 
data  reduction  is  presented. 

2.  4  CONSIDERATION  OF  THE  THRUSTOR  AS  A  PULSED  DEVICE 

In  the  previous  contract  period  it  was  shown  that  the  thrustor  can 
reach  quite  high  efficiencies  at  high  power  levels.  In  certain  applications 
the  total  developed  thruBt  is  not  necessary  but  the  high  efficiencies  are 
attractive.  This  leads  one  to  consider  the  possibility  of  using  the  thrustor 
in  a  pulsed  mode  where  full  thrust  is  developed  for  short  times;  i.  e.,  a 
duty  cycle  which  gives  the  desired  average  thrust,  but  with  an  efficiency 
commensurate  with  the  high  peak  thrust  values. 

We  consider  here  the  use  of  a  pulsed  mercury  vapor  feed  system 
and  we  choose  one  which  is  presently  undergoing  development  since 
pulsed  vapor  systems  are  not  easy  to  come  by  and  this  one  appears  to  be 
promising.  The  thrustor  has  never  been  operated  in  the  mode  investigated 
here  so  some  broad  assumptions  are  implicit  in  these  considerations, 
namely  that  we  may  pick  out  performances  in  the  steady- state  and  apply 
them  to  the  pulsed  mode.  The  particular  operating  conditions  are  those 
shown  in  Figure  B-12  of  the  Report.  The  pulsed  feed  system  is  that 
under  development  at  the  General  Electric  Space  Sciences  Laboratory4. 

The  device  (called  the  Exploded  Liquid  Metal  (ELM)  Injector)  produces 
a  pulsed  mercury  vapor  stream  and  has  been  tested  over  a  range  of  duty 
cycles  and  with  variable  mass/pulse.  Estimates  have  been  made  of  the 
actual  length  in  time  of  the  vapor  pulse  by  probing  downstream  of  the 
output  nozzle.  The  parameters  to  be  chosen  are 

a  -  mass /pulse 

4Gorowitz  B.  ,  Karras  T.  ,  Cook  C.  ,  and  Gloersen  P.  .  Summary  Report 
Research  and  Development  on  the  Pulsed  Plajma  Accelerator,  April  1967. 
Contract  NAS3-8906. 
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b  -  repetition  rate 
c  -  number  of  injectors 
d  -  pulse  length 

Since  wc  wish  to  use  the  data  of  Figure  B-1Z  we  must  choose  among 
the  parameters  n  a  reasonable  way  to  get  the  instantaneous  mass  flow 
rate  of  that  data.  Table  VII  of  the  GE  report  gives  us  the  following: 

.  9 

a  -  mass/pulse  =  l .  5X10  Kgm 
b  -  Repetition  Rate  =167  PPS 

The  thrustor  requires  a  good  degree  of  azimuthal  symmetry  in  the  mass 
flow  pattern  and  the  simplest  way  is  to  use  an  appropriate  number  of 
injectors  in  an  azimuthal  array.  A  reasonable  number  for  the  present 
thrustor  size  is  around  ten  injectors. 

c  -  number  of  injectors  =10 

The  pulse  length  is  estimated  in  the  report  as  "a  few  tenths  of  milliseconds" 
so 

d  -  pulse  length  =  0.  2  X  10  seconds. 

Combining  a,  c,  and  d  we  have 

m=  average  mass  flow  rate  during  the  pulse  on-time 
_  9 

m  =  f  ^ -■  ^  X  10  =  7.  5  X  10  **  Kgm  / sec 
.  2  X  10 

The  data  of  B-12  was  taken  at  7.  6  X  10  **  Kgm/sec.  The  pulse  on-time 
of  0.2  X  10"^  sec  is  about  50  cycles  of  the  RF  at  the  present  operating 
frequency  of  240  Kilohertz  and  this  seems  long  enough  to  assume  steady 
operation  of  the  thrustor.  The  instantaneous  mass  flow  rate  for  the  ELM 
will  diffei  fiOm  the  average  during  the  on-iime  and  actually  the  thrustor 
will  slide  through  its  performance  curve  from  higher  I  than  that  computed 
from  the  average  to  lower  but  since  the  time-shape  of  tne  mass  flow  rate 
is  not  known,  the  average  m  (averaged  only  over  the  on-time)  will  \  e  used. 
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The  thrustor  operating  point  is  chosen  to  be  at  an  efficiency  of  70% 
and  therefore,  from  Figure  B-12,  an  output  velocity  of  2.  5  X  10^  m/sec. 
The  instantaneous  power  input  is  33  KW  (23  KW  beam  power)  and  an 
instantaneous  thrust  of  1. 9  newton.  The  duty  cycle,  D,  of  the  ELM  at  a 
repetition  rate  of  167  PPS  is  .  033. 

The  average  thrust  <T>  is  given  by 

<T>  =  DT  where  T  is  the  instantaneous  thrust 
<T>  =  .  063  newton  (.  014  lbs) 

The  average  power  <P>  depends  on  whether  or  not  the  RF  is  pulsed. 
To  see  whether  it  is  feasible  to  leave  the  RF  on  continuously  (avoiding 
switching  circuits)  we  define  an  average  efficiency  <r|>  by: 


<n>  = 


_  u 

T 

2<rh>  <P> 


where  these  averages  are  taken  over  the  duty  cycle: 

<P>  =  PD  +  P„  (1-D) 

where  P  is  the  instantaneous  power,  and  P  is  the  coil  ohmic  dissipation. 
Writing 


PR  -  aP 


we  have 


<n>  = 


2m  P  D  + 


D 

a  (1-D) 


D  +  a  (1-D) 


where  r\  =  70%  at  the 
operating  point  and 


Therefore  <h>  =  “  23% 


This  appears  to  be  a  large  penalty  to  pay  in  order  to  avoid  pulsed 
P.F.  Using  pulsed  RF  the  average  efficiency  can  be  almost  as  high  as 
the  instantaneous  efficiency. 
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For  example  if  the  RF  is  pulsed  on  for  a  half  on-time  before  the 
vapor  pulse  and  left  on  for  a  half  cn-time  after  the  pulse  then: 

<ri>  =  1]  __1 — 

1  +  Q 

and  since  a«  1 

<*!>  =  i\  (1  -  a) 

In  this  case  then  <i"|>  =  65%  and  the  average  power  <P>  would  be  given  by 
<P>  =  PD  (1  +  a)  =  1.  2  KW 

The  average  thrust  can  be  lowered  by  lowering  the  repetition  rate  (the 
other  parameters  are  being  held  fixed  4t  the  known  values  which  work) 
as  can  the  average  power.  Since  the  repetition  rate  used  here  is  the 
highest  one  reported  it  would  not  be  realistic  to  discuss  raising  it. 

To  summarize  these  sample  calculations  we  have 

<T>  -  .  063  newton 

I  =  2500  seconds 
sp 

-6 

<m>  =  2.5X10  Kgm/sec 
<  P  >  =  1 . 2  KW 
<1  >  =  65% 
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3.  TKRUSTOR  VACUUM  CHAMBER 


3.  1  GENERAL 

The  major  changes  made  in  the  thruetor  chamber  system  was  the 
installation  of  a  KDH-1  30  forepump  and  a  completely  new  system  of 
water  baffles.  Otherwise,  the  chamber  is  as  described  in  the  Report. 

3.2  VACUUM  SYSTEM 

The  chamber  pumping  system  consists  of  a  number  of  freon-cooled 
baffles  for  pumping  mercury  vapor  and  an  MHG-900  mercury  diffusion 
pump  and  associated  forepumps.  During  the  past  contract  it  was  seen 
that  experiments  were  being  affected  because  the  diffusion  pump  became 
foreline-preesure  limited  during  the  experimental  runs.  At  the  begin¬ 
ning  of  this  contract  a  Kinney  KDH-1 30  pump  was  installed  in  parallel 
with  the  original  KP-15.  No  further  problem  of  foreline  limitations 
were  encountered.  In  addition,  the  increased  speed  substantially  de¬ 
creased  the  pumpdown  time.  Figure  3-1  is  a  chamber  pumpdown  curve. 

3.3  WATER  BAFFLES 

The  middle  water  baffle  system  shown  in  the  Report  {Figure  4-31) 
was  scrapped  and  a  new  system  built.  The  new  system  better  protects 
the  chamber  ports  from  direct  radiation  and  impingement  of  hot  gas 
(one  chamber  port  cracked  in  an  extended  run  during  the  previous  con¬ 
tract  period)  and  permits  a  calorimetric  power  measurement  over  the 
entire  chamber  surface.  The  baffles  are  shown  in  the  photograph  of 
Figure  3-2.  The  barrel  baffle  containes  four  iris  baffles  which  act  as 
radiation  and  impingement  barriers  and  the  five  sections  are  connected 
in  series.  The  chamber  water  systems  are  shown  in  Figure  3_3.  The 
No,  3  system  cools  the  upper  baffle  alone.  The  No.  2  system  supplies 
the  catcher,  lower  barrel  and  middle  barrel-iris  combination,  all  con¬ 
nected  in  series.  The  triangles  of  Figure  3-3  indicate  thermocouple 
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Figure  3-1.  Pumpdown  Curve  of  Chamber  with  KPH  130  Pump 


Figure  3-2.  New  Water  Baffles 


Figure  3-3.  Thermocouple  Locations  in  Water  System 


locations.  The  thermocouples  are  connected  to  measure  the  temperature 
drops  across  the  various  elements.  The  flow  in  the  water  system  is 
measured  by  diverting  the  flow  into  containers  where  the  time  to  fill  a 
fixed  volume  is  measured. 

The  pressure  in  the  system  was  continuously  monitored  by  an  ion 
gauge  connected  to  one  channel  of  the'Sanborn  recorder. 
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INSTRUMENTATION  AND  CALIBRATION 


* 


4.  1  GENERAL 

The  primary  parameters  that  describe  the  performance  of  the  accel¬ 
erator  are  thrust,  mass  flow,  and  electrical  input  power.  At  the  end  of 
Contract  AF33(657)-869’7  the  instrumentation  for  the  control  and  measure¬ 
ment  of  mass  flow  was  entirely  satisfactory,  but  that  for  the  measure¬ 
ment  of  thrust  and  electrica1 2  power  could  still  be  improved.  Accordingly, 
both  iivstems  were  redesigned  at  the  beginning  of  the  contract  period. 

The  thrust  platform  was  redesigned  on  the  basis  of  the  operational  ex¬ 
perience  with  the  previous  system.  A  new  power  measuring  system  was 
designed  on  the  basis  of  the  principles  developed  under  AFOSR  Contract 
AF49|638)-1251. 

4.  2  THRUST  STAND 

The  test  facility  for  the  thrustor  development  was  originally  designed 
to  accommodate  a  thrustor  firing  in  a  vertical  attitude,  with  the  thrust 
vector  pointing  up.  The  evolution  of  the  thrust  stand  retained  this  con¬ 
figuration  which  provides  two  important  advantages: 

1.  Effective  balancing  against  the  effect  of  mechanical  noise 
(vibration)  on  the  thrust  measurement. 

2.  Ease  and  precision  of  calibration. 

The  element  of  the  suspension  mechanism  are  shown  in  Figure  4-1, 

The  thrustor  is  attached  to  the  bottom  surface  of  the  platform  proper 
and  transmits  the  thrust  to  the  strain  gauge  cell  by  rnfano  Gf  aR  axial  steel 
rou  guided  by  the  two  ball  bushings.  The  static  weight  of  the  thrustor, 
platform  and  any  suspended  accessories  is  halanced  by  three  counter¬ 
weights  (one  is  shown  in  Figure  4-1  where  a  three -fold  symmetry  about 
the  axis  is  implied).  Exact  balance  can  be  obtained  by  adjusting  the 
position  of  the  fulcrums,  mounted  on  three  sliding  carriages  that  can 
be  simultaneously  moved  with  smoothness  and  precision  from  outside 
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of  the  vacuum  envelope.  It  is  clear  that  when  the  platform  is  statically 
balanced  it  is  also  balanced  against  vibration,  because  the  moments  about 
the  fulcrum(s)  are  equal  whether  they  are  due  to  the  acceleration  of  gra¬ 
vity  or  to  any  other  linear  acceleration  of  the  suspension,  except  for  the 
possible  elastic  flexing  of  the  structure  This  is  avoided  by  making  the 
structure  very  rigid  and  by  shock  mounting  the  entire  thrust  stand:  the 
resonant  frequency  of  the  total  mass  on  the  shock  mounts  is  distinctly 
different  and  much  lower  than  the  resonant  frequency  of  the  balanced 
mass  on  the  fulcrum(s),  so  that  the  vibration  of  the  vacuum  chamber, 
due  to  pumps,  refrigerators,  people  walking  about,  etc.  ,  is  attenuated 
to  an  entirely  acceptable  level  by  the  time  it  causes  any  forces  to  be 
applied  to  the  strain  gauge  cell.  Typically,  the  residual  detected  noise 
is  of  the  order  of  0.  003  newton  rme  when  measuring  a  thrustor  of  the 
order  of  0.  5  newton  with  a  suspended  thrustor  mass  of  5  Kg. 

The  strain  gauge  cell  (Daytonic  152-A)  was  found  to  give  very  accu¬ 
rate  and  reproducible  readings  when  operated  in  the  hostile  environ¬ 
ment  of  the  chamber,  which  normally  contains  mercury  vapor  at  a  pres¬ 
sure  of  the  order  of  10  ”4  torr.  The  ceil  is  in  direct  thermal  contact  with 
a  water  jacket  which  keeps  the  strain  gauge  cell  at  a  constant  temperature 
Water  cooling  of  strategic  elements  of  the  thrust  stand  prevents  the  gen¬ 
eration  of  drifts  of  the  indicated  zero  thrust  due  to  the  thermal  expansion 
of  the  suspension  and  balance  system  even  when  the  thrustor  components 
attached  to  the  platform  are  at  temperatures  in  excess  of  600°C. 

A  complex  problem  is  presented  by  the  flexible  leads  which  must 
carry  the  power  and  the  cooling  water  to  the  thrustor  without  generating 
any  spurious  forces  on  the  thrust  platform.  The  difficulty  of  the  problem 
can  be  appreciated  if  one  considers  the  various  requirements: 

1.  Current  capacity  to  100  ampere  rms  at  240  KHz 

2.  Voltage  insulation  to  5  KV  rms  when  immersed  in  a  rarefied 
conductive  hot  plasma. 

3.  Internal  pressure  to  80  psi  of  cooling  water 

4.  Indifference  to  thermal  expansion  and  local  electromagnetic 
forces 

The  basic  structure  of  the  flexible  leads  is  shown  in  Figure  4-2. 

The  current-carrying,  pres  sure -containing  elements  are  beryllium 
copper  bellows  which  provide  the  required  flexibility.  They  are  brazed 
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to  intermediate  copper  cylinders  which  are  perforated  in  a  hexagonal 
pattern  to  permit  the  passage  of  the  water  through  six  holes  and  the  pas¬ 
sage  of  a  thin  stainless  steel  cable  in  the  center.  The  axial  steel  cable 
absorbs  the  axial  tensile  stress  due  to  the  internal  pressure;  the  bellows 
are  stabilized  against  snake  instability  by  proper  choice  of  length/ diameter 
ratio  as  well  as  by  the  molded  silicone  rubber  jacket  which  also  supplies 
the  needed  electrical  insulations. 

Changes  of  temperature  of  the  cooling  water  cause  changes  of  length 
of  the  assembly  which  could  be  translated  into  spurious  thrust  if  the  re¬ 
sulting  forces  have  a  component  parallel  to  the  sensitive  axis  of  the  plat¬ 
form.  This  is  avoided  by  orienting  the  axis  of  the  leads  perpendicular  to 
the  axis  of  the  platform,  as  shown  in  Figure  4-3.  If  this  precaution  is  not 
observed,  the  indicated  thrust  is  subject  to  a  small  drift  with  a  time  con¬ 
stant  corresponding  to  the  thermal  constant  of  the  water-cooled  portion  of 
the  thrustor,  the  drift  becoming  noticeable  when  the  thrustor  is  fired 
periodically  on  and  off.  Figure  4-4  is  the  reproduction  of  two  traces  of 
the  thrust  channel  on  the  Sanborn  recorder  under  similar  firing  conditions 
(about  50  percent  duty  cycle,  one  minute  period,  3  to  4  kilowatt).  One 
trace  was  taken  with  the  leads  slightly  misadjusted  a  small  angle  from  the 
perpendicular,  the  other  with  the  leads  properly  adjusted. 

The  strain  gauge  is  calibrated  by  lifting  small  weights  which  normally 
rest  on  the  upper  surface  of  the  thrust  platform  (see  Figure  4-1).  The 
weights  are  in  the  form  of  concentric  rirgs  nesting  around  the  central 
vertical  rod  of  the  platform,  so  that  their  center  of  gravity  is  coincident 
with  the  axis  cf  the  rod.  They  are  lifted  by  remote  actuators  so  that  ex¬ 
actly  known  amounts  of  equivalent  thrust  can  be  introduced  at  any  time 
during  a  test  run,  whether  the  thrustor  is  firing  or  not.  The  weights  are 
machined  to  provide  exact  calibration  forces  of  0.  1,  0.  2,  0.  5,  1.  0,  2.  0 
and  5.  0  newtons;  by  lifting  two  together  0.  3,  0.  7,  3.  0  and  7.  0  newtons 
can  be  generated.  Calibration  over  the  entire  range  shows  that  the  re¬ 
sponse  of  the  strain  gauge  cell  is  extremely  linear  so  that  intermediate 
values  can  be  read  with  confidence.  The  maximum  probable  total  error 
in  any  thrust  reading  is  estimated  to  be  less  than  0.  007  newtons. 

The  thrust  platform,  once  balanced,  does  not  suffer  from  measurable 
zero  drift.  The  zero  drift  attributable  to  mechanica.1  vagaries  of  the  plat¬ 
form  is  less  than  0.01  newton  in  24  hours,  and  is  verifiable  at  any  time 
by  a  built-in,  remotely  actuated  locking  mechanism  which  leeks  the  plat¬ 
form  with  zero  force  applied  to  the  strain  gauge  cell.  The  drift  of  the 
baselines  as  recorded  on  the  Sanborn  trace  is  attributable  to  the  electronic 
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circuits  which  read  the  strain  gauge  output,  and  can  be  calibrated  out  by 
use  of  the  zero-force  platform  lock.  The  drift  rate  of  the  electronic  cir¬ 
cuits  is  of  the  order  of  1  newton/24  hours  and  is  somewhat  sensitive  to 
ambient  temperature.  However,  no  difficulty  results  in  obtaining  accurate 
readings,  because  the  operation  of  the  thrustor  is  automatically  interrupted 
for  two  seconds  every  minute,  thus  providing  a  current  zero  baseline  for 
all  recorded  quantities  . 

4.  3  POWER  MEASURING  SYSTEM 

The  electric  circuit  of  the  thrustor  is  a  simple  inductor  coupled  to  the 
plasma  circuit  as  shown  in  Figure  4-5,  where  the  equivalent  resistance  R 
represents  a  power  sink:  the  power  delivered  to  the  plasma  is  used  to 
ionize  the  gas,  heat  it,  and  accelerate  it.  The  effect  of  the  mutual  induc¬ 
tance  M  is  to  reduce  the  effective  inductance  of  the  thrustor  coil  as  well 
as  to  couple  in  a  resistive  component  R'  which  appears  in  series  with  the 
inductance. 

The  operation  of  the  thrustor  depends  on  the  variable  magnetic  field 
generated  by  the  coil,  and  is  therefore  strongly  dependent  on  the  current 
flowing  in  the  coil.  It  is  convenient  to  cancel  the  reactance  of  the  coil  by 
means  of  a  series  capacitance  so  that,  at  resonance,  an  adequately  heavy- 
current  can  be  made  to  flow  in  the  circuit.  The  capacitor  must  be  variable 
to  tune  the  circuit  as  the  effective  inductance  of  the  coil  changes  due  to 
the  plasma  loading  (as  a  function  of  mass  flow  and  power  level).  The 
extent  of  the  detuning  caused  by  the  plasma  is  reduced  by  adding  an  ex¬ 
ternal  inductor  in  series  which  the  circuit,  which  acts  as  ballast.  In 
addition,  a  small  inductor  is  added  for  sensing  the  current  flowing  in  the 
circuit.  The  complex  circuit  then  appears  as  in  Figure  4-6. 

Ll  and  R^  indicate  respectively  the  inductance  and  resistance  of  the 
current  sensing  inductor,  L,2  and  R£  those  of  the  thrustor,  L3  and  R3 
those  of  the  ballast  inductor,  while  L  and  R  represent  the  plasma  loading 
coupled  to  the  thrustor. 

Voltage  Vj  is  used  to  measure  the  current  Ijj 

v,  =  Ri  V  +  jwLi  rL 

Voltage  Vt?  is  measured  across  the  thrustor  and  contains  an  inductive 
component  juL^  II  plus  a  resistive  component,  representing  the  ohmic 
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losses  of  the  thrustor  coils  plus  the  power  delivered  to  the  plasma  which 
appears  as  an  equivalent  resistance  R': 

^  2  =  V1  +  ■*U)L2  +  *L  ^  + 

The  series  resonant  voltage  Vg  i'->  the  sum  of  all  inductive  and  resis¬ 
tive  voltages  minus  the  voltage  drop  across  the  tuning  capacitor,  which 
cancels  the  inductive  components,  as  indicated  by  the  vector  diagram  in 
Figure  4-7.  From  this  diagram  it  is  clear  that 

Vs  =  :L  <R1  +  *2  +  R'  *  R3> 

when  Cg  is  adjusted  so  that  all  inductive  components  are  exactly  cancelled. 
This  condition  is  identified  by  two  indications:  the  phase  angle  0  is  90°, 
and  the  absolute  magnitude-  of  Vg  is  a  minimum.  Then  the  power  dissipated 
in  the  circuit  is 


W  =  Vi.  *  :L2  <R1  +  R2  +  R1  +  R3>- 

Now  R}  is  very  small  and  negligible  in  comparison  with  the  other 
resistance,  R2  and  R1  are  the  quantities  of  interest  which  must  be  mea¬ 
sured,  and  R3  (the  resistance  of  the  ballast  coil)  must  be  known  and  sub¬ 
tracted  from  the  total  to  obtain  the  power  actually  delivered  to  the  thrustor. 

The  value  of  R3  at  240  KHz  was  accurately  measured  by  a  calorimetric 
method:  the  actual  power  dissipated  in  the  coil  was  measured  in  the  cooling 
water  for  a  series  of  values  of  1^  up  to  40  amperes,  and  the  equivalent  R3 
computed,  with  the  result 


R^  =  0.  281  ohrn. 

Therefore,  the  power  delivered  to  the  thrustor  is 

WS-  VSIL-0'“1IL2 

It  should  be  noted  that  the  power  dissipated  in  R3  does  not  appear  in 
the  measurement  of  power  with  the  It^V^  cos  0  method,  because  the  phase 
of  Vl  is  measured  at  the  terminal  of  L2  and  is  not  affected  by  R3. 
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The  most  sensitive  indicator  of  the  correct  tuning  of  Cg  was  found  to 
be  the  phase  angle  0.  Both  0  and  9  are  instrumented  and  recorded  on  two 
traces  of  a  Sanborn  recorder;  a  reading  is  considered  valid  only  when  0 
is  close  to  90° . 

4.4  ERROR  ANALYSIS 

Two  independent  methods  are  now  available  for  the  measurement  of 
the  power  delivered  to  the  load:  the  product  of  Vg  and  1^,  and  the  product 
of  VL,  IL,  and  cos  9-  The  current  IL  is  common  to  both  methods,  but 
the  other  quantities  are  separately  measured  by  independent  instruments. 
The  probable  accuracy  of  the  power  measurement  can  be  computed  as 
follows: 


For  the  series  resonant  method  we  have,  taking  into  account  the 
correction  for  the  resistance  of  the  ballast  inductor, 

ws  =  vsil-°'281Il2 

Accordingly,  the  absolute  value  of  the  expected  maximum  value  of  the 
error  in  Wg  is 

I  AWgl  =  Vg  |AlJ  +  II  I  AVgl  -  0.  562  IL  |  AlJ 

=  (Vg  -  0.  562  IL)  |AlJ  +  IL  |  AVgl 

The  electrical  power  is  also  measured  as  the  product  of  the  coil 
current,  the  coil  voltage  and  the  cosine  of  the  phase  angle 


=  VT  IT  cos  <p  =  V,  IT  sin  v|j 
9  L  L  L  L 


where  ^  =  90°  -  <p  is  the  angle  actually  measured  by  the  phase  meter  as 
it  is  inserted  in  the  circuit.  In  a  similar  manner,  the  absolute  value  of 
the  expected  maximum  error  in  W ^  is 

|  AW  |  =  V_  sin  ^  I  Al  |  +  I  sin  ^  \  AV  j 
9  S  L  L  L 

+  V  IT  cos  vp  !  A^| 

1_j  L 
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If  equal  weight  is  given  to  the  two  measurements  of  electric  power, 
the  most  probable  value  is  the  mean  of  the  two 


W 

e 


W  + 
s  <p 

2 


and  the  maximum  probable  error  is  the  rms  of  the  maximum  expected 
contributions  due  to  the  independently  measured  parameters:  Vg,  V-^, 
1^,  ip .  Note  that  Ij^  is  common  to  both  expressions,  whereas  the  other 
parameters  are  not.  Therefore  the  contribution  due  to  the  maximum 
expected  error  in  Ij^  is  the  mean  of  the  contributions  to  the  Wg  and  W^, 
cases  taken  before  squaring. 


w  =  — 

e  /I 


|l^  sin  I  AVj  +  ILVL  cos  ^  I  A'I'l]' 


+  iiL4vsr+Ksin'i'+vs-  0-562  v  |mji: 


1/2 


The  limits  of  measurement  accuracy  can  be  taken  as: 


|AVT  |  =  0.01  VT  +  20  volt 
L  L 

|  AVI  =  0.01  Vc  +  2  volt 

|  Al  |  =  0.  01  I  +  0.  2  ampere 
L  L 

|  AipJ  =  0.  5  degrees 

and  the  true  value  of  We  can  be  expected  to  be  between 

W  -  AW  <  W  <  W  +AW 
e  e  e  e  e 


:  AWC  is  an  upper  limit  to  the  error  and  can  be  taken  as  perhaps  three 

times  the  standard  deviation  of  the  measurements.  We  can  indicate  for 
brevity  the  upper  and  lower  limits  of  the  measurements  with 
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•-’•jv  ^-s .it&jSf*, 


wj  =  we  +  4  We 
w;  =  We  -  a  We 

and 


W"  <  W  <  W+ 
e  e  e 


Similarly  we  can  compute  the  upper  and  lower  limits  for  the  measurement 
of  the  mean  exhaust  velocity  v,  the  beam  kinetic  energy  and  the 
efficiency  rj: 


v  =  T/m 


V  <  V  <  v  * 


~±  — 
V  =  V 


Wk  =  T  /2m 


w,  <  w,  <  w, 

k  k  k 


W?  =  Wu 
k  k 


>*V(¥)W 


ti  =  T2/2fn  W 


*1  <*!  <  *1 


± 

*1  = 


AT  can  be  taken  to  be  0.  03  newton  and  Am/ m  is  essentially  negligible, 
however,  it  has  been  taken  as  0.  01. 


The  computation  of  the  v-or  and  lower  limit  of  each  performance 
parameter  for  each  set  of  measurement  involves  much  laborious  arithmetic. 
The  small  digital  computer  available  in  the  laboratory  has  been  programmed 
to  compute  and  print  out  the  follow  ing  quantities: 


-+ 
v  , 


w. .  w. .  w  .  w ..  w 


k’  V 


s' 


0’ 


*1 


given  the  inputs  of  the  measured  m,  T>  Vs>  VL>  ^  and  ijn  Thus,  it  is 
possible  to  obtain  the  upper  and  lower  limit  of  each  performance  parameter 
within  minutes  after  measuring  the  data. 

It  is  also  possible  to  verify  at  a  glance  that  Ws  and  W^,  two  independent 
measurements  of  the  same  quantity  We,  fall  within  the  expected  maximum 
error  of  each  other  and  of  their  mean  We. 

It  should  be  noted  that  the  instrumentation,  the  measurement  and  the 
computation  are  simple  and  straightforward  if  the  electrical  quantities  to 
be  measured  are  sine- waves  of  constant  frequency.  However,  the  plasma 
coupled  to  the  inductance  of  the  thrustor  is  a  nonlinear  load  and  harmonics 
are  generated  in  t1  -t  circuit.  The  series-resonant  circuit  presents  a  low 
impedance  to  the  fundamental  current  for  which  it  is  tuned,  but  a  high 
impedance  to  the  harmonic  frequencies;  therefore  even  a  small  percentage 
of  harmonic  current  in  the  coil  causes  an  appreciable  distortion  of  the 
voltage  waveform  appearing  at  the  Vg  terminals.  Furthermore,  a  distor¬ 
ted  waveform  creates  problems  ir  the  phase  measurements.  This  fact 
must  be  taken  in  account  and  a  suitable  filter  added  to  the  signal-processing 
circuits. 
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5. 


CONSTRUCTION 


5.  1  GENERAL 

The  copper  coil  must  be  protected  from  contact  with  the  mecury  vapor 
(plasma)  for  two  reasons: 

1.  The  mercury  vapor,  when  condensing  on  the  copper,  alloys  with 
it  resulting  in  quick  destruction  of  the  coil. 

2.  The  conductive  plasma  shunts  the  inductance  of  the  coil  and 
results  in  turn-to-turn  and  turn-to-ground  arcs  which  also 
can  destroy  the  coil. 

The  protective  material  must  be,  therefore,  a  good  electrical  insul¬ 
ator,  compatible  with  mercury  vapor,  capable  of  withstanding  very  high 
temperatures  where  it  comes  in  contact  with  the  plasma.  In  addition,  it 
must  be  capable  of  forming  on  the  copper  coil  without  leaving  voids  and 
without  cracking  either  during  the  cure  or  during  operation.  A  general 
requirement,  also,  is  that  the  thermal  expansion  coefficient  of  the  insula¬ 
ting  material  should  match  that  of  copper,  if  possible;  or  otherwise,  that 
the  insulating  material  should  be  strong  enough  to  cause  plastic  flow  of 
the  (annealed)  copper  without  cracking. 

The  required  characteristics,  taten  together,  appear  to  describe  a 
castable  ceramic  of  which  several  are  on  the  market.  A  preliminary 
screening  of  those  listed  in  the  commercial  literature  restricted  the 
choice  to  several  sold  by  Sauereisen  and  by  Aremco;  a  final  screening 
based  upon  tests  of  samples  eliminated  all  but  the  Aremco  Ceramacast 
505.  This  material,  although  superior  to  all  other  available  materials 
(including  the  PT  209  ceramic  used  in  Contract  Ai  33(657)-8697)  under 
the  conditions  of  the  tests,  still  was  not  entirely  satisfactory  under  actual 
service  conditions;  but  it  did  provide  a  starting  point  for  the  development 
of  a  material  composition  and  a  casting  process  that  eventually  produced 
satisfactory  coils. 
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The  firat  problem  was  shrinkage  during  setting  and  cure.  Shrinkage 
is  a  common  characteristic  of  all  physical  systems  undergoing  a  harden¬ 
ing  transformation.  When  a  rigid  object,  such  as  the  copper  coil,  is  em¬ 
bedded  in  the  shrinking  material,  stresses  develop  which  may  or  m a/  hot 
result  in  cracks  depending  as  the  curing  process  develops  sufi'cient 
strength  after  or  before  most  of  the  shrinkage  takes  place.  In  the  case 
of  the  ceramacast  505  the  shrinkage  (about  1  percent)  occurred  before 
adequate  strength  was  developed,  and  fine  hairline  cracks  remained  in 
the  ceramic  after  cure.  The  cracks  extended  to  the  copper,  and  caused 
arcs  when  the  thrustor  was  fired.  A  first  and  basic  improvement  in  occur¬ 
rence  of  initial  cracks  was  obtained  by  the  use  of  a  pre-cast  and  pre-cured 
base  for  the  coil,  in  the  form  of  a  thin  ceramic  wafer  over  which  the  main 
casting  was  poured.  This  backbone  provided  a  nonshrinking  support  for 
the  ceramic  during  the  cure  as  well  as  a  flawless  surface  on  the  side  facing 
the  plasma.  Yet  cracks  could  still  develop  on  the  opposite  face  of  the  coil. 

The  problem  of  shrinkage  was  minimized  by  adding  to  the  composition 
an  inert,  granular  filler  (tabular  alumina)  so  as  to  occupy  most  of  the 
volume  with  a  non-shrinking  component  and  to  restrict  the  domain  of  the 
shrinking  component  to  the  interstices  between  granules  where  it  could 
perform  a  binding  function  without  changing  the  distances  of  the  granules 
greatly.  If  the  mixture  could  be  so  proportioned  that  the  granules  would 
fill  exactly  all  the  remaining  space,  a  non-shrinking  composition  would 
result. 

It  was  found,  however,  that  before  this  condition  could  be  obtained, 
air  would  remain  trapped  in  the  interstices  between  granules.  Minute 
air  bubbles  embedded  in  the  solid  ceramic  would  appear  at  first  sight 
essentially  harmless,  and  so  they  are  if  the  dielectric  stress  applied  to 
the  ceramic  is  dc  or  at  a  low  frequency.  If,  however,  the  applied  voltage 
is  at  a  sufficiently  high  frequency,  the  air  in  the  bubble  can  ionize.  The 
concentrated  heat,  possibly  combined  w;‘‘h  wall  bombardment  by  the  ions 
oscillating  in  the  applied  alternating  field,  slowly  enlarges  the  cavity  of 
the  original  bubble  and  extends  it  in  some  direction  presumably  deter¬ 
mined  by  local  weaknesses  of  the  ceramic  as  well  as  by  the  direction  of 
the  applied  field.  The  result  is  a  "worm  hole"  that  eventually,  after 
many  hours  of  slow  progress,  reaches  either  a  copper  surface  or  the 
external  surface  of  the  ceramic;  the  result  of  this  is  a  catastrophic  arc 
after  an  extended  period  of  satisfactory  operation.  Such  "worm  holes" 
were  found  in  all  the  coils  that  failed,  after  many  hours  of  operation, 
with  an  applied  voltage,  at  the  time  of  failure,  no  higher  or  appreciably 
lower  than  the  voltage  they  had  successfully  carried  at  various  times 
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before.  Some  coils,  to  be  sure,  failed  by  direct  breakdown  at  a  weak 
spot:  but  this  happened  quickly,  while  voltage  was  being  increased,  the 
first  time  that  sufficient  voltage  for  breakdown  was  reached.  In  all  cases 
of  direct  breakdown  there  was  either  a  crack  or  insufficient  thickness  of 
insulation  due  to  eccentricity  of  the  coil  in  the  mold. 

The  incidence  of  minute  bubbles  was  1  educed  by  mounting  the  mould 
on  a  vibrator  and  by  shaking  the  mixture  while  it  was  being  poured.  This 
procedure  did  not,  however,  permit  the  use  of  a  mixture  containing  enough 
granular  alumina  to  eliminate  shrinkage  without  retaining  a  few  bubbles. 

It  was  found  that  the  addition  of  a  small  quantity  of  sodium  silicate  to  the 
mixture  allowed  the  use  of  a  much  higher  percentage  of  granular  alumina, 
to  a  point  where  shrinkage  was  essentially  cut  to  zero.  The  sodium  sili¬ 
cate  apparently  acts  as  a  lubricant  at  the  points  of  contact  of  the  hard 
granules,  and  permits  closer  packing  under  vibration;  it  also  acts  as 
supplementary  binding  agent  between  the  solid  particles  and  increases 
the  mechanical  strength  of  the  casting. 

The  glassy  characteristics  of  the  sodium  silicate  have  a  detrimental 
effect,  however,  because  the  drying  of  the  casting  is  delayed.  Castings 
containing  only  Ceramacast  505  and  tabular  alumina  appear  to  have  such 
a  combination  of  porosity  and  capillarity  that  they  dry  essentially  uniformly 
throughout  their  cross-section.  A  series  of  sample  ingots  cast  together 
and  broken  at  various  intervals,  showed  that  the  interior  was  becoming 
progressively  and  uniformly  dryer  as  time  went  by,  so  that  they  could  be 
cured  in  the  oven,  as  specified  by  the  manufacturer  of  the  Ceramacast 
505,  about  36  hours  after  casting.  Similar  samples  containing  sodium 
silicate  developed  a  hard,  impervious  crust  tnat  kept  the  interior  quite 
wet  and  made  a  cure  impossible  for  many  days.  Attempts  at  curing  the 
casting  before  it  was  sufficiently  dry  resulted  in  the  evaporation  of  the 
liquid  inside  the  crust  with  consequent  "puffing"  of  the  coil  like  leavened 
bread. 

It  was  possible  in  the  end  to  find  a  compromise  between  the  percentage 
of  sodium  silicate  in  the  composition  and  the  drying/curing  schedule  such 
that  adequate  coils  could  be  produced  in  a  reasonable  time.  The  ceramic 
composition  was  further  improved  by  using  two  mesh  sizes  of  tabular 
aiumina  (analogous  to  the  sand  and  the  gravel  in  concrete)  and  by  adding 
a  small  percentage  of  chopped  quartz  fibers  for  checking  crack  propaga¬ 
tion.  Since  each  coil  required  several  weeks  for  its  fabrication,  and 
many  days  of  testing  before  its  quality  could  be  fully  assessed,  any  useful 
feedback  of  information  to  the  manufacture  of  coils  had  a  built-in  delay 
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of  several  weeks.  At  the  end  of  the  program  both  the  ceramic  composition 
and  the  manufacturing  process  were  adequate  but  not  necessarily  optimum. 
The  coils  made  near  the  end  of  the  program  survived  testing  at  a  sufficiently 
high  voltage  and  for  a  sufficient  length  of  time  to  fulfill  some  of  the  require¬ 
ments  of  the  experimental  thrustor  program.  The  ultimate  limits  of  coil 
quality  and  thrustor  performance  have  not,  however,  been  reached,  and 
further  progress  can  be  expected  from  a  continued  study  of  the  properties 
of  ceramic  aggregates  of  the  type  so  far  used. 

5.2  RADIATION-COOLED  COILS 

The  Statement  of  Work  required  an  investigation  of  the  possibility  of 
operating  the  thrustor  in  a  radiation -cooled  mode,  without  water  circula¬ 
tion.  There  are  some  basic  differences  between  the  radiation-cooled 
thrustor  and  a  water-cooled  one.  A  water-cooled  coil  consists  of  a 
copper  tube,  in  which  the  water  flows,  covered  by  a  layer  of  insulating 
material.  The  outer  surface  of  the  insulating  material  is  exposed  to  the 
hot  plasm:',  and  is  heated  by  radiation  and  plasma  bombardment.  The 
electric  current  heats  the  copper  wire  in  proportion  to  its  resistance. 

Under  steady  state  conditions  the  temperature  of  each  element  will  be 
such  that  the  total  power  received  is  equal  to  the  total  power  lost  (power 
and  heat  are  here  used  interchangeably).  Power  sinks  are  the  flowing 
water  and  external  space.  The  inner  wall  of  the  copper  tube  can  be 
assumed  at  the  temperature  of  the  water.  The  outer  surface  of  the  coil 
will  take  the  temperature  which  will  radiate  (according  to  AaT^)  the  net 
amount  of  power  difference  between  the  power  delivered  by  the  plasma 
and  the  power  delivered  by  conduction  to  the  copper  and  hence  to  the 
water;  therefore  the  temperature  prefile  inside  the  coil  will  be  as  shown 
in  Figure  5-1.  A  temperature  gradient  is  established  in  the  insulating 
material  with  the  highest  temperature  at  the  surface.  The  case  of  the 
radiation  cooled  coil  is  entirely  different.  The  only  sink  is  external 
space;  therefore  the  outer  surface  must  assume  a  temperature  at  which 
it  radiates  the  sum  of  the  power  input  from  the  plasma  and  the  power 
input  from  the  copper;  all  other  things  being  equal,  therefore  the  sur¬ 
face  temperature  must  be  higher  than  in  the  water-cooled  case;  in  addi¬ 
tion,  a  temperature  gradient  must  be  established  so  that  the  heat  developed 
in  the  copper  by  I^R  losses  can  be  delivered  to  the  surface,  and  this  is 
in  the  opposite  direction  to  that  of  the  water-cooled  case.  Therefore  the 
temperature  of  the  copper  must  be  higher  than  the  temperature  of  the 
surface.  The  situation  is  illustrated  in  Figure  5—1.  It  is  clear  that  in 
the  radiation-cooled  case  temperatures  are  everywhere  much  higher  than 
in  the  water-cooled  case,  so  that  one  likely  mode  of  failure  is  the  melting 
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Figure  5-1.  Temperature  Diagrams  of  Water-  and  Radiation -Cooled  Coils 
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of  the  copper;  in  addition,  the  stresses  on  the  insulation  are  such  that  it 
may  tend  to  crack  because  the  outer  layers  are  stressed  in  tension  by  the 
hotter,  expanded  inner  layers.  For  this  reason  it  is  important  that  the 
material  used  has  high  mechanical  strength  and  high  thermal  conductivity. 
High  thermal  conductivity  can  be  obtained  by  the  substitution  of  granular 
beryllia  for  the  granular  alumina  in  the  aggregate.  This  was  not  done, 
however,  during  the  past  period  because  beryllia  is  much  more  expensive 
than  alumina.  It  was  decided  to  carry  out  all  experimental  development 
work  with  the  cheaper  alumina  until  a  composition  is  found  to  be  mechan¬ 
ically  satisfactory  (minimum  shrinkage,  no  cracks,  no  air  intrusions, 
etc.).  At  this  point  beryllia  can  be  substituted  volume  for  voh’me  in  the 
aggregate.  The  shift  to  water-cooled  coils  made  this  phase  of  the  pro¬ 
gram  unnecessary.  To  summarize,  the  radiation-cooled  coil  has  in 
common  with  the  water-cooled  coil  the  following  modes  of  failure: 

1.  Immediate  breakdown,  upon  first  application  of  voltage  in  the 
case  of  pre-existent  cracks  or  pinholes  exposing  the  copper 
to  the  plasma. 

•  2.  Delayed  breakdown  due  to  air  inclusions  which  grow  to 

"worm  holes.  " 

3.  Direct  voltage  breakdown  when  the  voltage  gradient  exceeds  the 
dielectric  strength  of  the  insulation.  This  happens  when  the 
voltage  is  increased  for  the  first  time  to  a  sufficient  level  in 
the  case  of  thin  or  partially  cracked  insulation. 

In  addition,  the  radiation-cooled  co;1  has  the  following  peculiar  modes  of 
failure: 

4.  Cracks  may  develop  due  to  the  expansion  of  the  copper  and 
inner  ceramic  layers  against  the  outer  layers  as  the  internal 
temperature  increases. 

5.  If  a  sufficiently  high  temperature  is  reached  without  cracks, 
the  copper  may  melt. 

Early  experiments  were  carried  out  in  a  bell  jar  to  determine  the 
operating  temperature  that  could  be  expected  in  a  radiation-cooled  coil. 

A  number  of  coils  were  tested  with  low  frequency  (60  cycles)  power, 
whereby  voltage  and  current  could  be  accurately  measured  by  standard 
instruments,  and  surface  temperatures  could  be  measured  by  thermo¬ 
couples  without  the  problems  attendant  to  the  induction  heating  of  the 
thermocouples  themselves  by  high  frequency  fields.  The  results  are 
shown  in  Figure  5-2;  also  shown  for  reference  is  the  slope  proportional 
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to  T^.  It  appears  that  the  heat  loss  from  the  coil  was  entirely  governed  by 
radiation,  since  the  slope  of  the  observed  points  is  closely  parallel  to  the 
slope  of  T^.  A  computation  of  the  radiated  power  per  unit  area  of  coil  sur¬ 
face  showed  that  the  composition  used  (straight  Ceramacast  505)  had  an 
emissivity  close  to  that  of  a  black  body  for  the  infrared  wavelengths  in¬ 
volved.  The  experiment  also  showed  that  the  operating  temperature  of  a 
coil  could  be  expected  to  be  in  excess  of  600°K,  and  probably  close  to 
800 'K.  All  the  coils  tested  developed  cracks  that  would  have  led  to  failure i' 
of  type  4  if  the  applied  voltage  had  been  high  enough;  this  was  taken  as  an 
indication  that  the  ceramic  insulation  lacked  the  mechanical  strength  to 
withstand  the  forces  exerted  by  the  expanding  copper  (similar  tests,  re¬ 
peated  later  with  an  improved  ceramic  composition,  resulted  in  melting 
the  copper  tube  in  its  ceramic  envelope  without  cracking  the  ceramic;  the 
composition  used  was  that  tabulated  in  Table  5-1  opposite  coils  39  and  40). 

The  same  series  of  tests  was  used  to  measure  the  resistivity  of  the 
copper  when  hot.  The  results  are  shown  in  Figure  5-3j  also  shown  is  the 
resistance  change  as  given  by  the  handbook.  This  information  was  used 
for  computing  the  ohmic  losses  of  the  coils  in  Section  2.2 

During  the  course  of  the  development  program,  the  gradually  improving 
materials  and  techniques  shifted  the  normal  mode  of  failure  from  type  1  to 
type  2  to  type  3.  Type  4  may  have  occurred  once  or  twice.  Type  5  was 
not  observed;  but  tests  on  copper  tubes  embedded  in  ceramic  samples  of 
more  advanced  composition  were  melted  by  the  passage  of  heavy  currents 
without  breaking  the  ceramic  as  mentioned  above. 

Every  coil  was  given  a  serial  number  from  the  time  fabrication  was 
started,  and  a  detailed  record  of  every  step  of  the  procedure  was  kept 
for  reference  against  performance.  Because  of  the  long  delay  in  infor¬ 
mation  on  feedback,  however,  the  benefits  of  experience  could  not  be 
applied  to  serially  successive  coils,  and  blocks  of  coils  were  produced 
that  upon  test  behaved  in  rather  similar  fashion.  This  had  the  disadvan¬ 
tage  that  progress  was  necessarily  slow,  but  on  the  other  hand  the  advan¬ 
tage  that  a  small  amount  of  statistics  could  be  collected  to  prevent  undue 
influence  of  freak  failures. 


It  was  thus  found  that  inner  coils  were  much  more  vulnerable  than 
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Thus  it  happened  that  several  outer  coils,  produced  early  in  the  program, 
remained  available  as  spares  to  the  end,  and  were  never  used;  and  soon 
production  of  outer  coils  was  stopped  altogether  and  all  efforts  were 
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TABLE  5-1 

COMPOSITION  OF  COILS 


Serial 

No. 

Type 

Cooling 

Composition 

Failure 

Mode 

Life 

1 

2 

3 

4 

5 

6 

Hr 

Min 

See 

13 

Inner 

Radiation 

100 

80 

22 

I 

5 

19 

Outer 

Radiation 

100 

80 

22 

* 

S 

33 

45 

20 

Outer 

Radiation 

100 

90 

22 

2 

70 

15 

15 

21 

Inner 

Radiation 

100 

100 

23 

1 

8 

33 

1 

22 

Inner 

Radiation 

100 

100 

22 

2 

7 

.  42 

58 

23 

Inner 

Radiation 

100 

100 

22 

1 

50 

47 

24 

Inner 

Radiation 

100 

100 

22 

2 

1 

13 

37 

25 

Outer 

Radiation 

Dented  —  Not  completed 

26 

Outer 

Radiation 

100 

iOS 

22 

*** 

27 

Outer 

Radiation 

Dented  -  Not  completed 

28 

Inner 

Radiation 

100 

105 

22 

29 

Outer 

Radiation 

100 

170 

35 

10 

30 

Inner 

Radiation 

100 

170 

35 

10 

2 

1 

30 

23 

31 

Inner 

Radiation 

100 

170 

35 

10 

2 

1 

30 

44 

32 

Inner 

Radiation 

100 

170 

35 

10 

2 

2 

46 

54 

33 

Inner 

Radiation 

100 

170 

35 

10 

3 

10 

36 

34 

Inner 

Water 

70 

B5 

50 

33 

2 

59 

34 

18 

35 

Outer 

Water 

70 

95 

50 

6 

2 

31 

** 

100 

48 

26 

36 

Outer 

Water 

70 

120 

60 

6 

4 

30 

*** 

37 

Inner 

Water 

70 

95 

50 

6 

2 

31 

3 

14 

49 

21 

38 

Inner 

Water 

70 

95 

50 

6 

4 

31 

3 

3 

46 

39 

Inner 

Water 

70 

120 

50 

8 

4 

31 

♦ 

53 

03 

55 

40 

Inner 

F.adiation 

70 

120 

60 

9 

4 

30 

41 

lnne  r 

Water 

70 

130 

65 

10 

4 

29 

3 

33 

51 

29 

Ingredients 

1  Ceramacast  505 

2  Tabular  Alumina  14-24  Mesh 

3  Tabular  Alumina  60-80  Mesh 

4  Sodium  Silicate  Solution  Robinson  1583 

5  Chopped  Quartz  Fibers,  1/4"  Long 

6  Distilled  Water 


*No  failure,  retired  because  of  mechanical  reasons. 
**No  failure. 

***Not  used.  Available  as  a  spare. 
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devoted  to  the  production  of  inner  coils.  Consequently  the  outer  coils  in 
service  at  the  end  of  the  program  reflect  materials  and  techniques  not  as 
refined  as  those  used  in  the  corresponding  inner  coils.  Table  5-1  is  a 
summary  of  the  materials  used  in  the  series  of  radiation-cooled  coils  that 
were  tested  with  significant  results  (excluding  early  coils  that  were  sub¬ 
jected  to  preliminary  tests  only  and  not  used  in  an  operational  thrustor). 

The  quantities  shown  as  ingredients  are  proportional  units  of  weight  (grams) 
except  the  quantity  of  sodium  silicate,  which  is  given  in  milliliters  of 
saturated  solution  in  water. 

The  life  times  shown  in  the  last  columns  are  for  reference  only.  The 
lifetime  is  significant  only  when  considered  in  conjunction  with  the  power 
level  at  which  the  coil  was  operated,  and  in  particular,  with  the  current 
(determining  l^R  losses  anu  therefore  temperature)  in  the  case  of  radiation- 
cooled  coils,  and  voltage  (determining  diselectric  stress)  in  the  case  of 
water-cooled  coils.  This  information  is  presented  in  Table  5-2,  which 
shows  the  chronological  sequence  of  the  use  of  the  coils  in  an  operational 
thrustor.  Tables  5-1  and  5-2  show  in  their  proper  places  the  water- 
cooled  coils  which  were  built  in  sequence  with  the  radiation -cooled  ones. 

5.  3  WATER-COOLED  COILS 

During  much  of  the  program  the  emphasis  of  the  development  effort 
was  directed  towards  obtaining  steady-state,  long  duration  operation  of 
the  radiation-cooled  thrustor,  with  the  specific  objective  of  100  hours 
life.  Before  this  objective  could  be  reached  the  effort  was  re-directed 
towards  obtaining  the  highest  possible  efficiency  rather  than  the  longest 
life  time;  and  since  efficiency  could  be  made  higher  by  using  higher  cur¬ 
rents  in  the  coils,  water-cooled  coils  were  called  for  rather  than 
radiation -cooled  ones. 


The  severe  operating  conditions  of  the  radiation-cooled  thrustor  were 
such  that  the  materials  were  subjected  to  an  equivalent  accelerated  life¬ 
testing  for  the  structures  in  the  water-cooled  mode,  so  that  a  vastly  im¬ 
proved  technology  for  building  water-cooled  thrustors  was  available  when 
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water-cooled  mode  at  mere  than  twice  the  voltage  used  in  the  radiation- 
cooled  mode  (and  more  than  twice  the  current);  the  outer  coil  did  not  fail, 
and  of  the  two  inner  coils  used  one  did  not  fail  and  the  other  failed  by 
direct  breakdown  (mode  3)  at  a  point  of  weakness  due  to  an  assembly  error ^ 
which  resulted  in  a  localized  lack  of  insulation  inside  one  of  the  supporting 
ceramic  legs;  the  failure  happened  at  the  highest  voltage  ever  used  in  this 
series  of  tests. 


TABLE  5-2 


CUMULATIVE  RECORD  OF  THRUSTOR  FIRING 


Inner  Coil 

Outer  Coil 

Maximum 

Stress 

Date 

On  Time 

Serial 

No. 

On  Time 

Serial 

Nc. 

Hr 

Min 

Sec 

Hr 

Min 

Sec 

Ampere 

Kilovolt 

10  April 

22 

0 

05 

55 

19 

0 

05 

55 

1 1  April 

1 

05 

52 

1 

05 

52 

12.  1 

0.800 

12  April 

4 

20 

22 

4 

20 

22 

13  April 

6 

18 

51 

6 

18 

51 

12.0 

0.750 

14  April 

7 

22 

58 

7 

22 

58 

16.0 

1.200 

4  May 

7 

29 

54 

7 

29 

54 

18.2 

1.000 

5  May 

7 

31 

13 

7 

31 

16 

19.  5 

1,  1 40 

7  May 

7 

42 

58* 

7 

42 

58 

17.  9 

1.000 

12  May 

23 

0 

04 

20 

7 

47 

18 

22.  5 

1.  350 

14  May 

0 

17 

59 

8 

11 

07 

22.  3 

1.  300 

15  May 

0 

50 

47* 

8 

33 

45** 

22.0 

1.250 

22  May 

24 

0 

26 

00 

20 

0 

26 

00 

18.5 

1.  140 

23  May 

1 

13 

37* 

1 

13 

37 

17.  0 

1.060 

1  June 

30 

1 

30 

23* 

2 

44 

00 

16.0 

1.000 

6  June 

31 

0 

09 

30 

2 

53 

30 

17. C 

1.070 

7  June 

1 

3(' 

4 

24 

14 

16.  5 

1.000 

9  June 

32 

1 

4-. 

19 

6 

13 

33 

18.  0 

1.  120 

10  June 

2 

4  3 

54* 

7 

15 

09 

lb.  3 

1.050 

12  June 

33 

0 

10 

36* 

7 

25 

45 

1 9  June 

0 

31 

51 

7 

58 

36 

18.  3 

l.  130 

20  June 

3 

49 

05* 

ii 

14 

50 

17.  2 

1.  100 

24  June 

34 

1 

42 

47 

12 

57 

37 

17.2 

1.100 

25  June 

14 

32 

21 

25 

47 

22 

16.  B 

1.000 

26  June 

29 

55 

15 

41 

10 

16 

17.  5 

1.000 

27  Jui.e 

48 

45 

00 

60 

00 

01 

15.  4 

0.950 

28  June 

59 

34 

14 

70 

15 

55* 

14.5 

0.85U 

5  July 

59 

34 

18* 

35 

0 

00 

04 

7  July 

4 

16 

07 

4 

16 

11 

8  July 

37 

14 

49 

21* 

14 

49 

25 

17.  5 

1.  100 

9  July 

33 

3 

46* 

14 

53 

A 

11  July 

39 

2 

01 

04 

16 

54 

16 

12  July 

% 

20 

11 

47 

35 

04 

59 

23.  3 

1.410 

1  3  July 

42 

39 

34 

57 

32 

46 

27.0 

1. 670 

14  July 

53 

03 

45** 

66 

56 

57 

30.0 

1.850 

26  July 

41 

5 

04 

01 

72 

00 

58 

29.0 

1. 770 

27  July 

16 

16 

01 

83 

12 

58 

34.  0 

1. 980 

.’8  July 

27 

28 

01 

94 

24 

58 

37.0 

2.  130 

•J  rt  .  .  1 .  . 

a.  / 

*1  *3 

J  J 

e  i 
* 

9A* 

/  - 

1 00 
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The  procedure  for  making  a  water-cooled  coil  of  the  latest  type  com¬ 
prises  the  steps  listed  below. 

1.  Annealed  copper  tube  (0.  125"  o.  d.  ,  0.  012"  wall)  is  bent  and  formed 
to  shape  in  an  aluminum  mold  (Figure  5-4).  The  strains  and  the 
work  hardening  due  to  the  bending  of  the  tube  are  removed  by 
annealing  again  the  tube  in  the  mold. 

2.  The  formed  tube  is  removed  from  the  mold  and  0.  250"  diameter 
feed-lines  are  brazed  to  the  terminals  of  the  coils.  The  assem¬ 
bly  is  chemically  cleaned,  leak  checked  by  helium  leak  detector, 
pressure  tested  at  80  psi,  and  cleaned  again.  Coils  at  this  stage 
are  shown  in  Figure  5—5. 

3.  The  coils  are  dipped  in  a  slurry  of  Ceramacast  505  with  added 
sodium  silicate  to  build  up  a  first,  closely  adherent  layer.  The 
layer  is  dried,  cured  and  carefully  inspected  for  flaws,  breaks, 
bubbles  or  lack  of  adherence. 

4.  Inner  coils  are  small,  compact  and  self-supporting;  individual 
turns  do  not  tend  to  sag  and  touch  ea:  h  other.  Outer  coils  are 
limber  so  that  turns  may  touch.  Outer  coils  are  dipped  again 
and  the  process  repeated  until  a  sufficient  thickness  of  insula¬ 
tion  is  built  up  to  insure  turn-to-turn  insulation  during  the 
casting  process. 

5.  A  0.  060"  -  0.  080"  thick  ceramic  wafer  is  cast  in  the  bottom  of 
the  casting  mold  (see  Figure  5-6;.  The  wafer  is  vibrated,  dried 
24  hours,  and  cured  24  hours  at  95“ C. 

6.  The  coil  is  suspended  and  centered  in  the  mold  above  the  wafer. 
Precise  positioning  is  required  to  avoid  shifts  of  location  that 
result  in  thin  ceramic  walls  and  consequent  danger  of  direct 
breakdown. 

7.  Three  ceramic  (Mullite)  tubes  are  positioned  so  that  their  ends 
will  be  embedded  in  the  casting.  The  tubes  serve  as  mechanical 
supports  of  the  coil  in  the  thrustor  and  are  located  in  standarized 
positions  so  that  all  coils  are  mechanically  interchangeable. 

8.  The  ceramic  aggregate  is  poured,  vibrated,  air-dried  36  hours, 
cured  at  50°C  for  36  hours  and  at  95“C  for  48  hours. 

9.  The  coil  is  pulled  from  the  mold.  The  feedlines  are  trimmed 
to  length  and  coated  with  Silastic  RTV  891.  Three  coats  are 
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applied  at  intervals  of  no  less  than  ten  hours  to  insure  complete 
cure  of  the  preceding  layer.  Swagelock  fittings  are  attached  to 
the  feedlines. 

10.  The  coil  is  baked  1?.  hours  at  150DC. 

11.  The  coil  is  weighed  and  stored  with  dessicant  in  a  sealed  container 
until  needed. 

Figure  5-7  shows  an  inner  coil,  an  outer  coil  and  the  nozzle  mounted 
on  an  assembly  fixture  for  fitting  the  connecting  leads.  Figure  5- 8  shows 
a  complete  thrustor  mounted  on  the  thrust  stand. 
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6.  THRUSTOR  PERFORMANCE 


6.  1  GENERAL 


The  Statement  of  Work  required  demonstration  of  thermal  equilibrium 
for  the  validity  of  a  performance  measurement,  mainly  for  the  purpose  of 
insuring  that  no  mercury  could  condense  in  the  feed  system  or  evaporate 
from  a  previously  condensed  state  so  as  to  cause  the  thrustor  to  operate 
with  an  instantaneous  mass  flow  different  from  that  measured  at  the  meter¬ 
ing  feed  pump.  This  requirement  was  observed  by  noting  that  the  thermal 
time  constant  for  heating  and  cooling  processes  in  the  thrustor  system  is 
of  the  order  of  15  minutes  (temperatures  are  monitored  by  thermocouples 
at  six  points  of  the  mercury  feed  system  and  two  points  of  the  platform), 
and  by  requiring  no  measurable  change  of  temperature  anywhere  in  the 
system  for  at  least  one  hour  before  the  time  of  an  official  observation. 


It  was  observed,  however,  that  as  long  as  the  temperatures  v/ere 
everywhere  higher  than  the  condensation  temperature  of  mercury  at  the 
local  mercury  pressure,  so  that  no  condensation  could  occur  anywhere  in 
the  mercury  feed  system  between  the  boiler  and  the  nozzle,  then  the  mea¬ 
sured  thrust  and  electrical  loading  did  assume  their  final  value  within 
seconds  of  a  change  of  either  mass  flow,  as  determined  by  the  metering 
feed  pump,  or  of  applied  electric  power;  no  changes  could  be  measured 
thereafter  for  as  long  as  mass  flow  and  power  were  kept  constant.  Simi¬ 
larly,  it  was  observed  that  the  measured  values  were  not  affected  by  the 
duty  cycle;  the  readings  were  identical  whether  the  thrustor  was  fired 
for  two  seconds  every  minute,  as  it  was  done  at  the  beginning  of  a  run  for 
outgassing  the  thrustor  and  the  chamber,  or  for  thirty  seconds  every 
minute,  as  it  was  done  after  several  hours,  when  the  outgassing  process 
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when  the  chamber  was  fully  outgassed  (firing  was  interrupted  for  two 
seconds  every  minute  to  guard  against  drifts  of  the  zero  baseline  of  the 
performance  recording  system;  thus  the  highest  normal  duty  cycle  was 
96.  ?%). 
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It  appears  therefore  that  the  practice  of  Jong  firing  periods  at  constant 
power  and  mass  flow  before  a  measurement  could  be  accepted  as  valid  was 
extremely  conservative.  Nevertheless  this  practice  was  observed  in  all 
measurements,  because  thf*  concurrent  goal  of  lifetime  testing  required 
that  the  thrustor  be  continuously  fired  for  many  consecutive  days,  so  that 
time  was  available  for  long  waits  between  changes  of  operating  conditions. 
Only  at  the  end  of  the  radiation-cooled  phase,  when  water-cooled  coils 
were  to  be  tested,  the  firing  time  required  for  acceptance  of  a  measure¬ 
ment  was  slightly  shortened,  with  the  concurrence  of  the  Air  Force  Pro¬ 
ject  Officer,  in  order  to  provide  a  larger  number  of  data  points  in  the 
available  time,  since  the  emphasis  was  shifted  from  long  life  time  to 
higher  performance. 

No  difference  in  performance  was  observed  between  radiation- cooled 
and  water-cooled  thrustors  over  the  range  of  operating  parameters  ac¬ 
cessible  to  both  modes.  In  fact,  all  measurements  were  extremely  con¬ 
sistent  and  repeatable  from  unit  to  unit  whenever  the  units  were  operated 
at  the  same  mass  flow  and  power  level. 

The  range  of  operating  parameters  accessible  to  the  thrustor  was 
explored  by  maintaining  a  constant  mass  flow  and  varying  the  coil  cur¬ 
rent  II,  or  by  keeping  the  current  constant  and  varying  the  mass  flow; 
the  series  resonant  capacitor  Cs  was  always  tuned  to  provide  an  indicated 
phase  angle  9  =  90°  and  a  minimum  value  of  V8.  In  general,  for  each 
new  coil,  input  power  was  increased  in  steps,  starting  at  a  low  value  {of 
the  order  of  two  or  three  kilowatt)  which  was  maintained  during  outgassing 
of  the  thrustor  and  chamber,  and  then  increasing  the  coil  current  in  steps 
of  two  or  three  amperes  :ms;  the  mass  flow  was  changed  in  steps  of  no 
more  than  10%  and  each  setting  was  maintained  through  at  least  one  tim¬ 
ing  cycle  of  the  automatic  mass  flow  measuring  apparatus  (see  Report, 
Section  5.  3.  1). 

This  procedure  resulted  in  a  detailed  map  of  the  performance  param¬ 
eters  over  a  range  of  mass  flow  between  0.  035  und  0.  16  grams/ sec  and 
power  levels  up  to  about  12  kilowatt.  It  was  found  that  the  most  sensitive 
indicator  of  correct  tuning  was  the  phase  angle  9  and,  to  a  slightly  smaller 
extent,  the  phase  angle  i|i.  Both  angles  responded  promptly  and  distinc¬ 
tively  to  small  changes  of  mass  flow.  Near  the  end  of  the  program,  after 
several  thousand  hours  of  operation,  the  metering  pump  in  the  mercury 
feed  system  showed  some  signs  of  wear:  the  cyclic  cam  lobe  which  drives 
the  pump  pistons  developed  a  flat  spot  where  the  cam  followers  had  some 
play.  This  play  resulted  in  a  small  irregularity  in  the  delivery  of  the 


mercury  to  the  boiler  coincident  with  the  passing  of  each  piston  from 
suction  to  delivery.  The  small  transient  periodic  irregularities  were 
detectable  in  the  recorded  traces  of  both  phase  angles  but  not  in  the  re¬ 
corded  trace  of  the  developed  thrust  nor  of  the  other  electrical  parameters 
(voltages  and  currents).  The  effect  was  particularly  obvious  at  the  very 
low  mass  flows,  where  the  transient  irregularities  represented  a  larger 
proportional  fraction  of  the  total  flow. 

All  performance  parameters  were  continuously  recorded  on  two 
8 -channel  Sanborn  recorders.  Recorder  No.  1  recorded  the  following 
channels: 

1.  Thrust  (newton) 

2.  Coil  current  II  (ampere  rms) 

3.  Coil  voltage  Vl  (kilovolt  rms) 

4.  Cos  0  (-1  <  cos  0  <  +1) 

5.  Series  resonant  voltage  V8  (kilovolt  rms) 

6.  Phase  angle  i|<  (degrees) 

7.  Plate  resonant  voltage  of  transmitter  Vb  (kilovolt  rms) 

3.  Plate  current  of  transmitter  Ig  (ampere  d.  c.  ) 

Recorder  No.  2  recorded  the  following  channels: 

1.  Chamber  pressure  by  ion  gauge  (10-3  Torr  full  scale) 

2.  Foreline  pressure  by  Thermocouple  gauge 

3.  Plasma  light  output  by  photomultiplier 

4.  Analog  signal  of  Cs  capacitance 

5.  Spare 

6.  Eight  multiplexed  channels  of  water  temperatures 

7.  Eight  multiplexed  channels  of  water  temperatures 

8.  Eight  multiplexed  channels  of  thrustor  temperatures 


107 


i 


The  normal  operating  procedure  during  the  tests  was  as  follows: 

Every  twelve  hours,  at  midnight  and  noon,  the  paper  was  cut  and  the 
roll  removed  from  each  Sanborn  recorder,  and  the  data  were  analyzed 
during  the  following  twelve -hour  period. 

The  chart  from  Recorder  No.  2  was  examined  first,  and  the  following 
data  transcribed  on  specially  prepared  forms  and  or.  standardized  graphs: 

1.  On-off  ratio  (duty  cycle)  when  less  than  96.  7%. 

2.  Total  cumulative  on  time  for  the  coil  set  assembled  in  the  thruster. 
These  were  taken  from  Channel  3,  which  indicated  whether  the 
plasma  was  actually  ionized  or  not. 

3.  Thrustor  temperatures. 

4.  Capacitance  of  tuning  capacitor  Cs. 

5.  Notes  and  comments  written  on  the  chart  by  the  operator  at  the 
time  of  occurrence. 

Inspection  of  the  graphs  showed  the  time  periods  in  which  the  tempera- 
tv.re3  and  other  conditions  were  constant,  for  the  required  long  intervals, 
and  such  times  were  noted  as  potential  sources  of  valid  performance 
points. 

The  chart  from  Recorder  No.  1  was  examined  next,  and  the  drift  of 
the  thrust  baseline  was  transcribed  on  the  same  graphs  that  were  show¬ 
ing  the  thrustor  tempeiatures  (it  was  thus  found  that  the  long  term  drift 
of  the  thrust  baseline  did  not  correlate  with  the  thrustor  temperature,  out 
with  the  ambient  tem-ierature  and  was  due  to  the  electronic  circuits  in 
the  Instrument  rack).  Then  the  time  periods  previously  noted  for  their 
unii-».m  condition  were  examined  in  close  detail,  and  if  the  electrical 
parameters  were  found  to  be  equally  uniform*  the  values  were  transcribed 
in  another  form  as  a  potentially  valid  performance  point.  The  final  accept¬ 
ance  of  the  readings  :ame  after  computation  of  the  e^ror  limits  and  veri¬ 
fication  of  the  power  readings  Ws  and  W^.  In  this  i  anner  anywhere 
between  one  and  eight  points  were  collected  each  day. 


6.  2  EXPERIMENTAL  DATA 

Table  6-1  shows  the  measured,  valid,  performance  points.  The 
points  are  arranged  in  order  of  increasing  mass  flow,  and  for  each  mass 
flow,  in  order  of  increasing  coil  current.  The  event  identification,  listed 
for  reference  purposes,  is  built  up  oy  the  month,  day,  (decimal  point), 
hours  and  minutes  of  the  record.  The  other  columns  are  labeled  by  the 
symbols  used  in  Section  4.  4. 

The  information  contained  in  Table  6-1  can  be  usefully  presented  in 
graphical  form  by  transcribing  the  data  in  a  suitable  format.  Figures  6-1, 
6-2,  and  6-3  present  the  value  of  as  a  function  of  T  for  three  approximate 
values  of  the  mass  flow  m.  Each  point  is  represented  by  a  rectangle 
limited  by  the  maximum  probable  errors  defined  in  Section  4.  4,  i.  e.  , 
the  top  and  bottom  margins  of  the  rectangle  are  q+  and  tj-,  and  the  right 
and  left  margins  of  (T  +  0.  01)  and  (T  -  0.  01)  newton. 

From  the  point  of  view  of  a  potential  user,  the  most  interesting  pre¬ 
sentation  is  that  in  the  (T,  v)  plane,  where  performance  parameters  are 
directly  related  to  the  desired  mission  requirements.  Since  v,  the 
average  exhavst  velocity,  is  defined  as  v  =  T/ m,  it  is  proportional  to 
the  specific  impulse  (ISp  =  v/g)  so  that  the  v  axis  can  also  be  labeled 
with  corresponding  values  of  Isp.  By  definition,  the  lines  of  constant  m 
are  fixed  in  the  (T,  v)  plane  and  are  shown  for  reference.  The  following 
illustrations  show  contours  of  constant  coil  current  1^,  coil  voltage  Vl, 
input  power  We,  and  efficiency  q.  When  plotting  the  experimental  points 
on  the  graphs  it  was  found  that  the  amount  of  smoothing  required  to  gen¬ 
erate  regular  curves  was  extremely  small.  In  practically  no  event  did 
a  parameter  fall  more  than  one  increment  width  away  from  the  ave’rage 
locus  of  the  parameter;  for  example,  all  points  with  II  =  35  ampere  fall 
between  the  locus  of  II  =  ampere  and  that  of  II  =  37  ampere. 

The  trends  which  are  apparent  in  the  illustrations  are  unmistake  able; 
but  the  area  so  far  covered  lies  below  the  transition  to  the  high  ionization 
efficiency  observed  at  high  power  level  under  Contract  AF33{657)-8697. 
Direct  extrapolation  of  the  data  collected  under  that  contract,  with  a 
co-rection  applied  for  the  lack  of  ferrites  in  the  current  configuration, 
indicated  that  a  change  of  ionization  efficiency  can  be  expected  in  the  range 
of  coils  current  50  <  Ij_,  <  90  ampere.  Currents  in  this  range  could  not 
be  reached  with  the  radiation-cooled  coils,  but  were  being  approached 
with  water  cooled  coils  at  the  end  of  the  contract  period-  Point  No.  03 
was  measured  with  the  highest  current  (43  amperes),  and  was  the  last 
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Figure  6-2.  Efficiency  vs  Thrust  for 


of  the  aeries;  it  was  also  obtained  with  the  lowest  mass  flow  (3.  5  X  10“-* 

Kb/ see).  Inspection  of  Figures  6-4  through  6-7  shows  that  this  point  does 
not  belong  to  an  orderly  extrapolation  of  the  other  measured  values,  but 
rather  indicates  some  change  of  the  functional  law  of  dependence  of  per¬ 
formance  on  current,  voltage  and/or  velocity.  Figure  6-8  shows  effici.  cy 
-r  as  a  function  of  current  1^  for  a  mass  flow  of  3.  5  and  3.  8  X  10“-*  Kg/  sec, 
with  the  data  points  represented  by  the  usual  rectangles.  The  peculiar 
position  of  point  No.  03,  above  the  orderly  extrapolation  of  the  other  mea¬ 
sured  values,  may  indicate  the  approach  of  a  transition  to  higher  efficiency; 
however,  in  view  of  the  results  of  Section  2.  3,  it  appears  more  likely  that, 
from  this  combination  of  mass  flow  and  current,  the  threshold  could  be 
reached  more  easily  by  an  increase  of  mass  flow  rather  than  an  increase 
of  current.  This  indeed  was  contemplated  in  the  plans  for  the  continued 
testing  of  the  coil,  but  was  prevented  by  the  failure  of  the  coil.  The  con¬ 
ditions  corresponding  to  point  No.  03  were  maintained  for  one  hour  and 
three  minutes;  the  inner  coil  failed  at  a  weak  spot  when  an  attempt  was 
made  to  increase  the  coil  current  prior  to  a  new  scan  of  the  mass  flow 
range  from  a  minimum  to  a  maximum. 

Figure  6-9  is  a  reproduction  of  the  traces  recorded  during  the  ob-  , 
servation  of  point  No.  03.  The  duty  cycle  was  temporarily  reduced  to 
50%  to  see  whether  a  better  vacuum  in  the  chamber  would  make  any  dif¬ 
ference  (it  did  not;  all  readings  taken  below  5  X  10“^  torr  were  indepen¬ 
dent  of  pressure).  For  comparison,  also  shown  are  the  traces  of  point 
No.  70  in  which  a  similar  value  of  thrust  was  obtained  with  a  much  larger 
mass  flow  at  a  much  lower  efficiency. 

Figure  6-9  also  illustrates  the  effect  of  the  mass  flow  fluctuations, 
due  to  the  play  in  the  feed  pump,  on  the  electrical  quantities.  Since  the 
mass  flow  is  the  lowest  ever  used,  and  the  wear  in  the  pump  is  the  worst 
(being  this  point,  03,  the  latest  in  the  series),  the  induced  fluctuations  in 
current  and  voltage,  as  well  as  phase,  are  quite  noticeable.  Capacitor 
C8  was  adjusted  to  cu.ae  out  the  series  resonant  circuit  under  conditions 
of  average  mass  flow,  as  indicated  by  0  =  90°  (cos  9  =  0).  This  condi¬ 
tion  happens  to  be  closely  obtained  during  the  firing  period  indicated  by 
an  asterisk,  and  the  reported  readings  are  taken  there.  It  is  of  interest 
to  note  that  the  phase  angle  9  is  a  sensitive  indicator  of  mass  flow  load¬ 
ing;  a  swing  in  one  direction  indicates  a  mass  flow  lower  than  average, 
in  the  other  higher,  as  shown  on  the  trace  of  cos  9.  Inspection  shows 
that  when  the  mass  flow  drops  below  average,  current  and  voltage  climb 
above  average,  and  the  thrust  record  shows  fast  fluctuations  as  if  condi¬ 
tions  in  the  thrustor  were  at  the  threshold  of  an  instability.  This  is  not 
the  case  for  point  70. 
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Figure  6-4.  Contours  of  Constant  Current  1^  in  the  (T,  v)  Pla.ne 
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2.  4  kilovolt  rma 


Figure  6-5.  Contours  of  Constant  Voltage  V.  in  the  (T,  v)  Plane 
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Figure  b-b.  contours  of  Constant  Power  Vvre  in  ilie  (T,  v)  Fletne 
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Figure  6-9.  Recorded  Traces  of  Observed  Performance 
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6.  3  MASS  UTILIZATION 


The  loss  of  efficiency  from  the  theoretical  values  can  he  attributed 
either  to  excess  power  required  to  ionize  the  mercury  vapor  above 
10,  4  electron-volts  per  ion  pair,  or  to  lack  of  ionization  for  part  ot  the 
effluent  mass,  or  to  a  combination  of  both  mechanisms. 

If  we  assume  that  the  total  mass  flow  m  is  divided  in  two  components, 
the  ionized  fraction  m  and  the  neutral  fraction  if^,  we  can  define  the  mass 
utilization  r|m 


m.  m.  . 

l  l 

”1  =  ~r~  -  - r — i” 

m  m  m.  +  m 
1  o 

Then,  if  we  assume  that  m^  is  accelerated  and  supplies  the  entire  mea¬ 
sured  thrust,  while  m,-,  is  not  accelerated  at  all  and  supplies  no  measur¬ 
able  thrust*,  we  can  compute  the  average  velocity  of  the  accelerated 
component 


v  =  T/m.  =  T/n  m,  =  v/q 
#  i  mi  m 


and  the  efficiency  of  the  acceleration  process 


q  =  T2/2  m.  W  =  TZ/2ti  m.  W  =  q/q 
v  i  e  'm  i  e  m 


If  the  mass  efficiency  can  be  independently  measured,  it  becomes  pos¬ 
sible  to  obtain  an  indication  of  the  relative  importance  of  the  two  aspects 
of  the  ionization  process;  how  many  ions  are  made,  and  how  much  power 
is  spent  in  making  them.  An  estimate  of  qm  is  obtained  by  plotting  the 
observed  points  on  a  graph,  derived  from  the  basic  theory,  in  which  the 
presentation  is  very  sensitive  to  rim  and  nearly  independent  of  ine  other 


*In  other  words,  if  we  assume  that  the  mean-free-path  X  for  charge 
exchange  collisions  is  large  compared  with  the  dimensions  of  the  inter¬ 
action  region.  The  information  in  Figure  2-45  shows  that  this  assump¬ 
tion  is  not  far  from  the  truth. 
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parameters.  Such  graph  is  found  in  the  form  of  a  family  of  curves  of 
constant  t)m  in  the  (T/Ii,  m)  plane,  already  shown  as  Figure  2-47, 

The  result  is  shown  as  Figure  6-10.  The  experimental  points  clus¬ 
ter  is  a  well  defined  band,  with  very  low  scatter,  with  approximate/y 
0.  05  <  qm  <  0.  1.  V/hereas  Figure  6-10  ;  rovides  an  immediate  visual 
presentation  of  the  distribution  of  rim,  Taole  6-2  shows  the  correspond¬ 
ing  value  for  each  observed  point,  identified  by  the  same  number  as  in 
Table  6-1.  Also  shown  are  the  values  of  v^  and  since  both  v^  and  t)* 
are  inversely  proportion?!  to  qm,  it  is  easy  to  compute  the  error  limits 
vj  and  v#  as  well  as  qj  and  qj  taking  the  non- starred  values  from  Table 
6-1  if  desired.  This  was  not  done  in  Table  6-2  for  simplicity).  In  addi¬ 
tion,  Table  6-2  shows  the  energy  per  ion  pair,  in  electron  volts,  corres¬ 
ponding  to  that  fraction  of  the  input  power  which  is  not  accounted  for 
elsewhere.  The  information  contained  in  Table  6-2  can  be  presented  in 
graphical  form  in  the  (T,  v*)  plane,  corresponding  to  Figures  6-4  through 
6-7,  except  that  now  the  velocity  axis  refers  to  the  velocity  of  the  ionized 
fraction  only  of  the  total  mass  flow  (see  Figures  6-11  through  6-14);  the 
presentation  of  Figure  6-11  is  of  particular  interest  because  it  affords  a 
direct  comparison  with  the  theoretical  predictions  of  Figure  2-44.  It  can 
be  seen  chat  the  agreement  is  quite  good.  By  inspection  of  these  figures 
one  comes  to  the  conclusion  that  the  thrustor  accelerates  very  effectively 
the  ioniaed  fraction  of  the  effluent  mass,  and  that  the  loss  of  efficiency  is 
predominantly  associated  with  the  escape  of  non-ionized  mass.  Also, 
the  ionization  mechanism  appears  to  be  effective  in  the  sense  that  not 
many  electron-volts  are  required  to  ionize  those  atoms  which  are  ionized. 
It  appears  therefore  that  a  substantial  increase  of  performances  can  be 
expected  by  the  action  of  any  mechanism  which  can  increase  the  number 
of  ions  present  in  the  interaction  region. 


The  existence  and  effectiveness  of  a  mechanism  of  this  sort  has  been 
observed  in  electrodeless  inductive  discharges  before;  not  only  under 
Contract  AF33(657)-8697,  but  also  during  a  project  sponsored  by  the 
Air  Force  Office  of  Scientific  Research  ^AF  49(638)-1251,  which  was 
concerned  primarily  with  making  accurate  electrical  measurements  of 

the  operating  power  and  the  resistance  ard  inductance  reflected  into  the 

J _ :  .  _i _  ..  1  _ i  t _ /  *  _  _i _ 5  _ s  j _ :  _ - •  .  .  _ _ 1  i _ At.  _ 
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5Penfold,  A,  S.  ,  and  Warder,  T,  C.  ,  "Research  on  Electromagnetic 
Plasma  Acceleration,  Volume  If,  An  Investigation  of  the  Various  Plasma 
Mischarges  Surrounding  a  Solenoidal  Coil  Excited  with  Current  at 
4  Megacycles,  "Final  report,  Contract  AF  49(638)-1251 ,  Litton  Pub¬ 
lication  5865,  Janu" ry  1  967. 
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TABLE  6-2 

PARAMETERS  OF  IONIZED  FLOW 


Event 

No. 

t/il 

•  * 

V 

c  -V 

ion 

01 

0.00575 

0.092 

64.62 

84.2 

781 

02 

0.00560 

0.087 

74.85 

88.8 

706 

03 

0.00744 

0.121 

74.87 

105.9 

04 

0.00656 

0.097 

48.70 

89.8 

269 

05 

0.00611 

0.091 

57.67 

87.5 

473 

06 

0.00621 

0.087 

62.67 

92.4 

324 

07 

0.00657 

0.093 

63.77 

95.8 

176 

08 

0.00620 

0.085 

48.12 

82.4 

496 

09 

0.00652 

0.083 

56.62 

89.5 

374 

10 

0.00799 

0.108 

49.21 

102.0 

11 

0.00719 

0.092 

42 . 62 

97.7 

42 

12 

0.00666 

0.087 

42.57 

88.3 

239 

13 

0.00740 

0.094 

46.44 

104.9 

14 

0.00616 

0.075 

49.50 

81.3 

565 

15 

0.00637 

0.077 

50.90 

83.3 

S19 

16 

0.00707 

0.084 

54.91 

94.4 

179 

17 

0.00722 

0.088 

54.69 

94.9 

160 

18 

0.00689 

0.074 

49.36 

92.2 

204 

19 

0.00702 

0.077 

61.45 

96.4 

.140 

20 

0.00754 

0.080 

45.08 

101.8 

21 

0.00786 

0.088 

43.02 

98.8 

21 

22 

0.00777 

0.085 

44.54 

101. 

23 

0.00750 

0.081 

27.35 

62.9 

442 

24 

0.00719 

0.074 

42.94 

103.5 

25 

0.00701 

0.073 

55.13 

91.6 

278 

26 

0.00662 

0.065 

48.70 

86.9 

358 

27 

0.00689 

0.068 

49.24 

89.5 

284 

28 

0.00791 

Q.  080 

45.24 

95.5 

95 

29 

0 , 00745 

0.071 

53.19 

96.5 

101 

30 

G. 00727 

0.068 

27.28 

53.3 

533 

31 

0.00624 

0.056 

27.34 

58.3 

535 

32 

0.00772 

0.071 

38.71 

93.4 

105 

33 

0.00765 

0.067 

30.32 

79.8 

233 

34 

0.00769 

0.069 

52.05 

94.9 

145 

35 

0.00757 

0.067 

53.61 

90.7 

294 

TABLE  6-2  (cont) 


Event 

No. 

T/IL 

V 

V* 

e-v 

ion 

_ 

0.00739 

0.064 

76.78 

118.5 

!37j  ■' 

0.00882 

0.076 

27.30 

83.6 

146 

38  * 

0.00789 

0.066 

31.25 

47.3 

091 

39 

0.00761 

0.064 

34.36 

49.9 

186 

40 

0.00857 

0.073 

28.20 

86.3 

126 

41 

0.00872 

0.072 

28.09 

85.5 

133 

42 

0.00818 

0.066 

36.38 

90.2 

144 

43 

0.00777 

0.064 

79.65 

132.5 

44 

0.00705 

0.055 

28.60 

70.2 

347 

45 

0.00750 

0,060 

25.92 

64.4 

373 

,  ■  46 

0.00764 

0.061 

54.98 

91.0 

298 

47 

0.00683 

0.053 

26.66 

74.6 

243 

48 

0.00773 

0.060 

27.62 

68.5 

351 

49 

0.00818 

0.064 

35.51 

53.7 

1090 

50 

0.00823 

0.065 

27.10 

70.6 

305 

51 

0.00823 

0.062 

27.95 

63.3 

453 

52 

0.00892 

0.067 

26.61 

92.3 

58 

53 

0.00731 

0.0S2 

27.24 

81.1 

173 

54 

0.00814 

0.059 

35.97 

91.5 

120 

55 

0.00854 

0.062 

56.66 

103.7 

56 

0.00795 

0.057 

52.92 

91.0 

276 

57 

0.00812 

0.056 

26.92 

90.7 

74 

58 

0.00818 

0.056 

55.61 

95.1 

157 

59 

0.00872 

0.063 

27.46 

90.3 

80 

60 

0.00877 

0.062 

27.80 

92.2 

65 

61 

0.01137 

0.081 

26.58 

115.1 

62 

0.01176 

0.085 

26.17 

121.3 

63 

0.01000 

0.077 

15.75 

37.5 

415 

64 

0.00859 

0.058 

20.64 

71.7 

168 

65 

0.01000 

0.069 

26.43 

102,7 

66 

0.00837 

0.054 

60.88 

101.2 

67 

0.00826 

0.054 

56.12 

96.2 

121 

68 

0.00914 

0.061 

27.32 

96.5 

26 

69 

0.00787 

0.050 

53.72 

86.7 

444 

70 

0.00914 

0.054 

55.67 

101.1 

m 


TABLE  6-2  (cont) 


Event 

No. 

T/Il 

V 

V* 

*1* 

e-v 

ion 

71 

0.00865 

0.051 

53.30 

92.9 

216 

72 

0.01172 

0.066 

25.34 

117.4 

73 

0.01233 

0.069 

24.07 

126.6 

74 

0.00845 

0.045 

67.79 

93.2 

244 

75 

0.00873 

0.047 

53.42 

93.5 

199 

76 

0.00799 

0.042 

24.07 

51.0 

558 

77 

0.00812 

0.043 

25.65 

53.9 

565 

78 

0.00874 

0.047 

25.40 

57 4* 

480 

79 

0.00416 

0.020 

21.17 

41.7 

628 

80 

0.01250 

0.072 

17.64 

83.4 

62 

81 

0.00799 

0.041 

32.34 

50.4 

031 

82 

0.00755 

0.038 

36.44 

42.7 

786 

83 

0.00965 

0.047 

23.61 

94.8 

30 

84 

0.01071 

0.055 

21.80 

97.0 

14 

85 

0.00918 

0.042 

29.89 

55.1 

729 

86 

0.00923 

0.042 

32.01 

63.6 

587 

87 

0.01098 

0.051 

29.30 

68.5 

396 

88 

0.01000 

0.038 

11.34 

62.3 

77 

89 

0.01117 

0.042 

27.76 

71.4 

310 
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Figure  6-11.  Contours  of  Constant  1^  in  the  (T,  v  )  Plane 
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Figure  6-13.  Contours  of  Constant  We  in  the  (T,  v  )  Plane 
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Army  Engineering  laboratory,  Fort  Belvoir,  Virginia  (DA-44-00c-AMC- 
1447  T),  in  which  the  emphasis  was  directed  toward  making  extensive 
diagnostic  measurement  of  the  plasma  state  using  spectroscopic,  micro- 
wave,  and  probe  techniques8.  In  all  cases  the  transition  was  observed 
to  be  abrupt  and  to  be  a  function  of  the  applied  drive:  when  the  applied 
current  and/or  voltage  and/or  magnetic  field  induction  exceeded  a  defi¬ 
nite  threshold,  then  the  gas  conductivity  increased  by  a  very  large  factor. 

The  available  evidence  seems  to  indicate  that  this  mechanism  can  be 
expected  to  begin  to  operate,  in  the  case  of  the  thrustors,  at  coil  currents 
above  approximately  70  amperes.  Since  coil  currents  of  this  magnitude 
have  not  been  reached  in  the  course  of  the  contract,  the  exact  location  of 
the  threshold  boundary  has  not  been  identified. 


8  Thornton,  J.  A. ,  "An  Investigation  of  the  Energy  Exchange  Processes 
in  Nonequilibrium  Plasmas,  "  Final  report,  Contract  DA-44-009-AMC- 
1447  (T),  Litton  Publication  5794,  January  1967. 
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